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Abstract 
We consider the Randall-Sundrum (RS) set-up to be a theory of flavour, as an alternative to Froggatt- 
Nielsen (FN) models instead of as a solution to the hierarchy problem. The RS framework is modified by 
taking the low energy brane to be at the GUT scale. This also alleviates constraints from flavour physics. 
^Z Fermion masses and mixing angles are fit at the GUT scale. The ranges of the bulk mass parameters are 

^H determined using a y^ fit taking in to consideration the variation in 0(1) parameters. In the hadronic 

O , sector, the heavy top quark requires large bulk mass parameters localising the right handed top quark close 

'^N to the IR brane. Two cases of neutrino masses are considered (a) Planck scale lepton number violation and 

LJ (b) Dirac neutrino masses. Contrary to the case of weak scale RS models, both these cases give reasonable 

^_^ fits to the data, with the Planck scale lepton number violation fitting slightly better compared to the Dirac 

case. In the Supersymmetric version, the fits are not significantly different except for the variation in tan/3. 
If the Higgs superfields and the SUSY breaking spurion are localized on the same brane then the structure 
of the sfermion masses are determined by the profiles of the zero modes of the hypcrmultiplcts in the bulk. 
^^ Trilinear terms have the same structure as the Yukawa matrices. The resultant squark spectrum is around 

^^^J ^2 — 3 TeV required by the light Higgs mass to be around 125 GeV and to satisfy the flavour violating 

OO constraints. 
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I. INTRODUCTION 

One of the celebrated solutions of the fermion flavour problem is the Proggatt-Nielsen Mechanism 
[1]. According to this prescription, the symmetry group of the Standard Model (SM) is augmented 
by a horizontal U{l)x group under which all the SM fermions and the Higgs field are charged. The 
effective theory includes a flavon field X and the Yukawa couplings are generated from the higher 
dimensional operators which are invariant under the U{l)x and the Standard Model (SM) gauge 
group. For example, the up-type quark mass matrix has the form: Yj^lj^j^^^^ '^"^ '^""QiH^Uj, 
where c/ is U{l)x charge of the / field and i,j are the generation indices. The flavon field X 
develops a vacuum expectation value (vev) such that 0.22 ~ Ac ~< X > /Mpi, Ac being the 
Cabibbo angle. Y-^ are taken to be 0(1) parameters. Fermion mass matrices including their 
mixing patterns can be fit to the data by choosing appropriate U{l)x charges for various fields. 
An UV completion of the model can be constructed by including heavy chiral fermions in to the 
theory; integrating these heavy fields would lead to the relevant non-renormalizable operators (for 
a review, see [2] ). 

The U{l)x symmetry introduces additional anomalies in to the theory and subsequently, strong 
constraints on the U{l)x charges for various fields. In supersymmetric models with a single flavon 
fleld, one typically has to resort to Green-Schwarz (GS) mechanism to cancel the anomalies. The 
solution set of U{l)x charges for the fermions and the Higgs which satisfy the fermion data as well 
as the anomaly cancellation^ requirement have been studied in [3-14] and recently updated in [15]. 
These models typically lead to large flavour violations at the weak scale in gravity mediated super- 
symmetry breaking models due to contributions from the U{l)x D-terms. While the constraints 
from the flavour sector on the available solutions are very tight, it may still be possible to ease 
them without requiring the superpartner masses to be very high [16, 17]. The flavour constraints 
may also be alleviated to some extent by considering U{1) x U'{1) class of models [18, 19]. In 
the present work, we will study the extra- dimensional alternative [20] to understand the flavour 
hierarchy in particular concentrating on the supersymmteric Randall- Sundrum (RS) set up. 

The Randall-Sundrum framework [21] which elegantly provides a solution to the hierarchy 
problem via warping in the extra dimensional space can also thought to be a theory of flavour. 
It has been observed sometime ago that the flavour changing neutral currents (FCNC) can be 
suppressed due to the so-called RS-GIM mechanism [22]. However, in the absence of additional 



^ However, with two singlet fiavons there exist a unique solution which is completely non anomalous [3] 



flavour symmetries the constraints from FCNC are still very strong( [23-27]) (Detailed analysis for 
the hadronic sector can be found in [28, 29] and references there in. For a recent thorough analysis 
in the leptonic sector, please see [30]). Given these strong constraints on the RS set up at the weak 
scale, one can ask the question whether RS is suitable to be a theory of flavour as well as a solution 
to the hierarchy problem simultaneously. It might be that RS as a theory of a flavour might be 
better suited at the GUT scale rather than at the weak scale. The Froggatt-Nielsen models are 
typically defined at scales closer to the Planck scale, so perhaps flavour physics might have its 
origins at the Planck scale. 

With this rationale, in the present work we will consider RS to span between the Planck and 
the GUT scales. The fermion masses are flt in terms of the bulk mass parameters of the various 
fields, which take the role of the U{l)x charges of the FN mechanism. However, these parameters 
are less constrained compared to the U{l)x charges, as no additional conditions such as anomaly 
cancellations are required on them. While this has been the common understanding, in [31] it was 
pointed out that imposing unification conditions on gauge couplings in a theory with localization 
of fermions or hierarchical wave functions leads to strong constraints which are exactly in the same 
way as the Green-Schwarz anomaly cancellation conditions [32] ^ In the present setup we do not 
impose these conditions. 

Extra dimensions at GUT scale were considered in [33] while the RS version was considered 
by the authors in [34, 35] and later by the authors in [31]. Our work, however, is very closely 
related to the work of [36] who have done a thorough analysis of fermion mass spectrum, weak 
scale supersymmetric spectrum and fiavour phenomenology, assuming a particular Grand Unified 
Theory (GUT) model in such a RS setting. However, differences exist. In the present work we 
have not assumed any specific GUT model. Furthermore, we have used a frequentist approach 
to do the fermion mass fitting. While this makes it hard to directly compare the results between 
the two works, we hope they provide a complementary set of results. We also have taken in to 
consideration the constraints from neutrino masses and mixing angles which can have a significant 
effect on the lepton fiavour violation and slepton decays. 

The equivalent description for the RS set up in four dimensions can be thought of as a composite 
Higgs coupled to fermions with couplings which parameterise the 'partial compositeness' of the 
fermions[37]. In SUSY case, this partial compositeness can also affect the structure of the soft 



^ Typically applied in FN models, the Green-Schwarz anomaly cancellation conditions requires the anomaly factors 
to be in a particular ratio such that they are cancelled in String theory. 



masses. 

In the first part of our work, our aim lias been to provide a range of bulk mass parameters which 
fit the fermion masses and mixing patterns at the GUT scale. We believe this can be useful for 
model builders and other phenomenologists working in flavour physics and looking for an alternative 
to FN models. We have considered both supersymmetric as well as non supersymmetric versions 
of the RS framework at the GUT scale while fitting the data. The supersymmetric case has the 
added advantage that it could lead to observable signatures at the weak scale. We consider the 
case where SUSY breaking is considered to be on the same brane as where the Higgs is localized, 
which is the GUT brane. In this case, the sfermion mass matrices are determined by the zero mode 
profiles of the corresponding N=l superfields and thus the information of the fermion masses is 
propagated in to the soft sector. It is far more striking for the A-terms which follow the same 
structure as the Yukawa couplings. The spectrum is highly non-universal at the high scale, but, 
its pattern is constrained due to the ranges of bulk mass parameters which are in turn are fixed by 
their fits to fermion masses. The running effects make the diagonal terms large at the weak scale. 

The rest of the paper is organized as follows. In section II, we detail the RS setup we consider 
and derive the structure of the fermion masses. In section III we present the fermion mass fits 
and present the ranges for the bulk mass parameters for both the non-supersymmetric and the 
supersymmetric cases. In section IV we address the issue of supersymmetric breaking and derive 
supersymmetric spectrum for a particular supersymmetric breaking case. We end with summary 
and outlook in the last section. In Appendices A , Band C, we have presented plots relevant for 
fermion mass fits. 

II. RS AS A THEORY OF FLAVOUR 

The Randall-Sundrum frame work consists of two branes separated by an single warped extra 
dimension [21]. The line element for the RS background is given as 

ds^ = e-'^'^yri^^dx^dx'' - dy^ (1) 

where < y < ttR. Here y = is identified as the position of the UV brane and y = irR is the 
position of the IR brane. The scale associated with physics on the UV brane is Mpi while that on 
the IR brane is TeV. The solution to the hierarchy problem is achieved by exponential warping of 
scales i.e, Mpi = e^^"" M^cak, where kR ~ 0(11). 

In the modified set up we consider here, the scale associated with the IR brane is Mqut- This 



can be achieved by choosing kR ~ 1.5. We define the hierarchy between the scales, e, for this 
scenario to be 



10~^ (2) 



Mpianck 

We consider both supersymmetric and non-supersymmetric matter fields to propagate in the bulk. 
Localisation of the respective zero modes is dependent on the corresponding bulk masses. In both 
the cases, we assume that the Higgs field (two Higgs fields in the case of supersymmetric models) are 
localised on the GUT brane. We now proceed to briefly review the derivation of the mass matrices 
and their dependence on the zero mode profiles in both supersymmetric and non-supersymmetric 
cases. A couple of points are important to note at this juncture. Firstly, the typical lowest KK 
mass for a warped background is given as rriKK = e~^^'^k. In the present set up, the 'large' 
warp factor ensures the lowest KK modes are very heavy i.e, rnxK = eA^ ~ Mqut and thus are 
decoupled from low energy phenomenology. We do not consider their effects in this work for low 
energy phenomenology. Secondly, it turns out that the dependence of the zero mode mass matrices 
on the profiles is very similar in supersymmetric and non-supersymmetric cases. However, the 
fermion mass data at the high scale in the supersymmetric case could be different from the SM 
one, due to the dependence on tan/3 as well as the different RGE for the Yukawa coupling as we 
will discuss in the next section. 

A. Standard Model case 

The non-supersymmetric case or the Standard Model case has been studied in many works 
[23, 28-30]. The main difference in the present case is that while those studies have considered a 
kR ~ 0(11), while in the present case it is 0(1). We thus assume a grand desert from the weak 
scale to the GUT scale, where RS framework sets in. No attempt is made to solve the hierarchy 
problem, but the fiavour problem has a solution in terms of the localisation of the fields in an extra 
dimension at the GUT scale. We follow the notation of [30] and present the final formulae for 
the Yukawa mass matrices. The details of the KK expansion and the corresponding ortho-normal 



relations can be found in [301 and references therein. The five dimensional action has the form: 



'-' — '-'kill + OYuk -\- Ju ~^ ^higgs 

Skin = j d'^x j dy ^/^ {L{ilJ)-mL)L + E{ilt)-mE)E + ... ) 
5Yuk = I d'^x I dy 7=^ ( YuQUH + YdQDH + YeLEH^ 5{y - ttR) 
Sy = f d'^x f dy ^/^(-^LHLH or YnLNh) 5{y - ttR) (3) 

where we used the standard notation with the Q, U, D standing for the quark doublets, up-type and 
down type singles respectively, L and E, N stand for leptonic doublets and charged and neutral 
singlets respectively. H stands for the Higgs doublet with H = iaH*. We have suppressed the 
Higgs action in the above. Two specific ways for generating non-zero neutrino masses are considered 
(a) by a higher dimensional term localised at the GUT brane and (b) Dirac neutrino mass terms 
similar to the other fermions. A'^^ is the five dimensional reduced Planck scale ~ 2 x 10^® GeV. 

After the Kaluza-Klein (KK) reduction and imposing the orthonormal conditions, we can derive 
the 4D mass matrices for the zero-modes of the fermion fields. They have the form: 

where F stands for all the Yukawa matrices F = U,D,E and N, if the neutrinos have Dirac masses. 
Cj and c'- represent the bulk masses of the respective matter fields (second line of eq.(3)); defined 
as, for example, m^;- = C£;.A:. i,j denote the generation indices. If the neutrinos have Dirac masses 
then their mass matrix is given by Eq.(4). In case the neutrinos attain their masses through higher 
dimensional operator, the mass matrix is given by 

In Eqs. (4, 5), we have defined Y' = kY and k' = 2kK. These are dimensionless 0(1) parameters 
of anarchical nature^. Eq.(4) are used to fit all the fermion mass data at the GUT scale i.e, 
up and down type quark masses and the (Cabibbo-Kobayashi-Masakawa) CKM matrix, charged 
lepton masses, neutrino mass differences and the corresponding PMNS mixing matrix. In the case 
neutrinos get their masses through higher dimensional operator, Eq. (5) is used instead to fit their 
mass differences and mixing angles. 



Note that the Yukawa couphngs in Eq.(9) are dimensionful, with mass dimensions -1. 



B. Supersymmetric case 

In the supersymmetric case, the matter fermions are represented by hyper-multiplets^ propa- 
gating in the bulk. In terms of the 4D, N=l SUSY language, they can be expressed as two N=l 
chiral multiplets, $, ^'^. Following [38, 40], we write the 5D action in terms of two chiral fields with 
a (supersymmetric) bulk mass term to be 



85= f d^x I (fee-'^^y ($f <5 + ^^^"■f) + / Soe-'^^y^^ (dy + m$ - h^ $ 



(6) 



where M$ = c^k is the bulk mass. In writing the above, the radion field is suppressed by taking 
its vacuum expectation value, < Re[T) >= R. The super field <I>^ is taken to be odd under Z2. 
Thus, only <I> has a zero mode. Since we have a theory at the GUT scale, the KK modes can be 
considered to be decoupled from theory. In the effective theory, the profile of the zero mode of the 
<I> is determined by [38] 

5. - (^ - c) k^ /(°) = (7) 

Thus f^^' = e'a"'^' y . The superscript '"^ stands for the zero mode, which we will drop subse- 
quently^. In this effective theory, where the higher KK modes are completely decoupled, we can 
write the effective 4-D Kahler terms for the Z2 even zero modes as [31, 34, 35] 

where we have substituted for the profile solutions of (7). After integrating over the extra dimension 
y, the terms in Eq.(8) pick up a factor Zp = tt^ — \u {^^^~^ — l) where F = Q,U, D, L,E, 



(l-2cF)fe 



as 



before. We choose to work in a basis in which the Kahler terms are canonically normalized. We 
thus re-define the fields as $ — t- —r^=^. 

The effective four dimensional MSSM Yukawa couplings are determined from the superpotential 
terms written on the boundary. For Higgs localized on the IR brane, keeping only the zero modes 
of the chiral superfields, the effective four dimensional superpotential is given as [31, 38] 

+ e^l-^^^^'^ye^l-'^^^'yv^HoLiEj + . . ) 6iy - ttR) (9) 



'^ N=l Supersymmetry in 5D has the particle content of N=2 Supersymmetry in 4D. The hypers can be expressed in 

N=l, 4D language as two chiral superfields, where as the Vectors can be expressed as a vector and chiral superfield 

[38, 39]. 
"" In the component form, the scalar component and the fermion components of the chiral super field $ have different 

bulk masses. However, the solution for the profile for the scalar and the fermion components turns out to be the 

same. 



The Higgs fields are canonically normalized as H^^d ~^ ^^^^H^^d- In the canonical basis, after the 
fields have been redefined the fermion mass matrices can be derived from Eq. (9) to be 

{Mf\, = !^y/.e(i-^-^i)'=-C(c.)e(4) (10) 

where Ci,c', denote the bulk mass parameters for various fields, ^(cj) are defined in Eq.(4) The 
mass matrix in Eq.(lO) can be approximated as {MF)ij '^ ^7^C(l)e ~'^'^'^J' '"'^ where the ^(c) is 
absorbed into the 0(1) parameters Y' and is now collectively referred to as 0(1). This is true only 
as long as the c parameter lies between and 1. But as we have seen earlier, in some realistic cases 
especially related to neutrino masses and the top quark, the values of \c\ could be large to fit the 
data. Redefining the 0(1) Yukawa by absorbing the c parameters would shift the ranges of the c 
parameters far away from what they are, especially in the case where \c\ > 1. As in the SM case, 
we define dimensionless 0(1) Yukawa couplings as Y' = 2kY and are of anarchical nature. While 
the Dirac masses for the neutrinos have the same structure as the other fermion mass matrices, 
the higher dimensional operator has a different form determined by the super-potential term 



w(4) 



I'dydiy - 7ri?)e-3-fe) (^e^^l-^^^^^ye^l-'^J^^y ^HuHuLiLj') 



(11) 



The neutrino mass matrix in this case is given as 



where k' = 2kK is the dimensionless 0(1) parameters. The fermion mass matrices carry the 
same form as in the SM and the supersymmetric cases and thus their dependence on q and 0(1) 
parameters is the same. 



III. FERMION MASS FITS 

From the previous section, we have seen that in addition to the bulk mass parameters, the 0(1) 
Yukawa parameters also play a role in fixing the fermion masses and mixing angles. We fit the 
masses and the mixing angles of the quark sector and the neutrino sector at the GUT scale for 
both the SM and the supersymmetric cases. We will use a frequentist approach, i.e, we minimise 
the x^ function, which is defined as follows: 

^^ii f'-f' ^j (13) 



eory ^^ ^j^^ theory number for the j^^ observable and y^^^ 
quoted by experiments with a measurement uncertainty of aj. In the present case the theory 



where, y- ^"^^ is the theory number for the j^^ observable and y^^^ is its corresponding number 



parameters are just the bulk mass parameters and the 0(1) Yukawa entries in supersymmetric 
and non-supersymmetric cases. We define < x^ < 10 to be a good fit and we try to find regions 
in the parameters space of bulk mass parameters and 0(1) Yukawa parameters which satisfy this 
condition^. The 0(1) Yukawa parameters are varied between -4 and 4, with a lower bound of 0.08 
on \Y\ to avoid unnaturally small Yukawa parameters. As far as the bulk mass parameters are 
concerned, since they are given as ck, we prefer to vary the c parameters between —1 to 1, so as 
not to go beyond the 5D cut-off, k. This will remove any possible inconsistencies in the theory 
due to non-perturbative Yukawa couplings. However, as we will see it is not always possible to fit 
the data within this range of c parameters. We will mention the range chosen specifically for each 
case. 

The minimization of the x^ function was performed using MINUIT [41]. We can minimise 
the hadronic and the leptonic sectors independently as they are dependent on different sets of 
parameters which are uncorrelated. The methodology is similar to the ones used in fermion mass 
fitting in GUT models [42, 43] and also the one used in [30]. 

A. Standard Model (SM) Case 

In this section, we present the fits in the SM case. For the GUT scale values of the quark and 
lepton masses and CKM mixing matrices, we use the results of [44]. In the analysis of [44], two 
loop RGE have been used to run the Yukawa couplings of the up- type quarks, down- type quarks 
and charged leptons from the weak scale all the way up to the GUT scale. For the neutrino data 
we used the publicly available package REAP [45] to compute the high scale values. The masses of 
the SM fermions at the GUT scale used in our fits are presented in Table I. The CKM and PMNS 
mixing matrices are presented in Table II. 

1. SM Quark Sector fits 

The up and down mass matrices are given in terms of fermion mass matrix of Eq.(4). The 
theory parameters which are varied simultaneously to minimise the x^ ^^ Eq.(13) include: three 
CQ^, each of c„^ and c^^ and 18 0(1) Yukawa parameters. We would expect that the light quarks 



We will mention the results with lower x^ a-t relevant places. 



TABLE I: GUT scale masses of fcrmions for the SM case 



Mass 

(MeV ) 


Mass 
(GeV) 


Mass 
MeV 


Mass squared Differences 


TO„ = 0.48+0;?° 

ma = 1.14tl!i^ 

m, - 22+_l 


TO, = 0.235t°:°^^ 

TO, = i.ot°:Sl 
TOt = 74.01;^:? 


TO - n 4606+"°"°"°°°'' 
irie — u.^ioyo^Q gQQQooo4 

TO - qq i4+"°"°°°* 
in^ — ^^-^^-0.0000089 

rur = leSS.SStoig 


Amj^ = 1.5+°;2; X 10-4 
Am|3 = 4.6+°;^!^ X 10-3 



TABLE II: Mixing angles for the hadronic and the leptonic sector for the SM case 



mixing angles (CKM) 



^12 - U./;^b_o 00087 

Q — n nzLi c;+o.oooi9 
013 = 0.0035t°;l]°} 



Mixing angles (PMNS) 



0.59: 



0.015 



023 = 0.79t^;i2 



'13 



= 0.154 



+0.016 
-0.016 



would be localised close to the UV brane ( c > 1/2 ) and the heavy quarks close to the IR brane 
(c < 1/2). However, for this particular range of 0(1) Yukawa parameters, it is difficult to fit the 
data for |c| within unity. We thus enlarged the range for the c parameters. 

The range chosen for the scan of the c parameters chosen is: —2 < cq-^^q.^ < 4, —3 < cqj < 1 for 
the doublets. —2 < Cd^^d243 < 3-^' fo'^ ^^^ down type singlets and —2 < Cu^^u2 < 4, —4 < 0^3 < 1 
for the up type singlets. We fit the quark masses and the CKM mixing angles at the GUT scale. 
The top quark is definitely lighter at the GUT scale, but still we see that most of the points that 
fit the data lie outside of |c| < 1. This is evident from the the negative values of the cqj and Cf/3 
that fit the data. 

The regions of c parameter space which satisfy the constraint of < x^ < 10 for the chosen 
scanning range are shown in Fig.(l) in Appendix A and the ranges are outlined in Table [III]. We 
see that the first two generation bulk mass parameters are concentrated on the positive c values 
where as the third generation, the doublet and more so the right handed top is localised close to 
the GUT scale brane. Comparing these results with that of the normal RS, we find that the masses 
for the light quark fields can be fit with c ^ 0.6 — 0.7. This can be attributed to the large warping 
where 0.5 < c < 1 is sufficient to reproduce the masses for light quarks [46]. 
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TABLE III: Allowed range of c parameters in the SM case. These parameters satisfy < x^ < 10 for the 
SM case. The corresponding figure is 1 in Appendix A. 



parameter 


range 


parameter 


range 


parameter 


range 


CQ3 


[0,3.0] 
[-1.95,2.36] 

[-3,1] 


CD2 


[0.78,4] 
[0.39,3.02] 
[0.39,2.21] 


CC/2 
CC/3 


[-0.97,3.98] 

[-1.99,2.43] 

[-4,1.0] 



2. SM Leptonic Mass fits 



Unlike the quark case, the fits in the leptonic sector are far more difficult and more constraining 
due to the small mass differences and the large mixing in the neutrino sector. As mentioned, we 
will consider two different cases of neutrino masses while fitting the leptonic data, 
(a) LLHH higher dimensional operator 

Planck scale lepton number violation is an interesting idea which manifests itself with higher 
dimensional operator suppressed by the Planck scale. In four dimensions such an operator generates 
too small neutrino mases. It is typically used as a perturbation over an existing neutrino mass 
model [47]. If not, it needs an enhancement of 0(10'^ — 10^) to be consistent with the data. In 
the standard RS framework close to the weak scale with bulk fermions, this higher dimensional 
operator is still constrained however for different reasons. While the neutrino masses can be fit by 
placing the doublet fields L close to the UV brane, the charged lepton masses become very tiny 
unless the singlet fields (E) are placed deep in the IR[30]. This leads to inconsistencies in the theory 
with large non-perturbative Yukawa couplings. The question arises whether the situation repeats 
itself when we consider the modified RS setup. This can be checked as follows. The neutrino 
masses are generated by the higher dimensional operator as given in Eq.(3). The corresponding 
neutrino mass matrix is given by Eq.(5) while the mass matrix for the charged leptons is given by 
Eq.(4). 



For simplicity assume cl- = c^V i. For cl < 0.5 the mass matrix in Eq.(5) becomes 



nil, = K 



2eA 



:i - 2cl) 



(14) 



It is clear that cl ~ —4 is required to get neutrino masses 0(0.04) eV for a warp factor for e ~ 10~^. 
As Cl increases, beyond 0.5, this formula is no longer valid, the neutrino masses become smaller 
and hence do not fit the neutrino mass data with 0(1) Yukawa couplings. Thus a mildly negative 
Cl should be able to fit the data without large inconsistencies. 

11 



A second enhancement can also come from the k', which is the corresponding 0(1) Yukawa. 
With this in mind, we enhance the range of the scanning of the Yukawa couphngs from 0.08 to 4 
to 0.08 to 10. This would help us to accommodate cl values close to ~ —1. The final scanning 
ranges we have chosen are: the doublets (c^-) are varied between -1.5 and 0.5, while the charged 
singlets were scanned between and 4. The region of c values which give a good fit to leptonic 
masses, i.e, satisfying the constraint < x^ < 10, is presented in Table[IV]. The plots for these 
ranges of c values are presented in Figs. (2) in Appendix A. 

TABLE IV: Ranges for scanned regions of the bulk leptonic parameters for the LLHH in the SM case which 
satisfy < x^ < 10. 



parameter 


range 


parameter 


range 


Cii 


[-1.5,-1.15] 


ce^ 


[2.8,4.0] 


CL2 


[-1.5,-0.97] 


CE2 


[1.8,2.4] 


CL3 


[-1.5,-1.22] 


CEs 


[1.2,1.69] 



(b)Dirac type Neutrinos 

The case of Dirac neutrinos is interesting possibility though it requires imposition of a global lepton 
number conservation^. The running of the neutrino masses from the weak scale to high scale is 
different in this case. However with a normal hierarchy of neutrinos and low tan/3 the differences 
are insignificant [49, 50]. 

Assuming that there is not much of a difference for normal hierarchy, we choose the following 
scanning range for the c parameters. The doublets {cl-) and charged lepton singlets (c^;-) are 
scanned within the range -1 to 4.5, while the neutrino singlets were scanned in the range 3.5 
to 9. Such a large value of the bulk mass parameters for the singlets is needed to suppress the 
corresponding neutrino masses sufficiently. The 0(1) Yukawa parameters were varied between 
0.08 and 4. Comparing the results of Dirac neutrino mass fits with that of the weak scale RS 
models, [30], we find that the cn are roughly a factor 7 — 8 larger compared to the cjsr at the weak 
scale. This is purely because of the weaker warp factor we are considering in the present case. 
Increasing the range of the 0(1) Yukawa parameters would only make things worse. The ranges 
for the c values corresponding to SM fits with Dirac neutrinos case are presented Table [V]. The 



''' In fact, it is possible to hide lepton number violation in this case through a careful location of the right handed 
fermion fields [48]. We will not consider this case here. 
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plots for the c parameters are presented case in Fig[3] in Appendix A. 



TABLE V: Ranges for the scanned regions of the bulk leptonic parameters for the Dirac case which satisfy 
< x^ < 10 for the SM case. 



parameter 


range 


parameter 


range 


parameter 


range 


CLi 


[ -1,2.9] 


CE^ 


[0.39,3.62] 


CATi 


[5.29,8.97] 


CL2 


[-0.99,2.7] 


CE2 


[-1.0,2.63] 


CN2 


[5.31,8.99] 


CL3 


[-0.99,1.98] 


CEs 


[-0.99,1.93] 


CNs 


[5.12,8.97] 



B. Supersymmetric Case 

The analysis for the case with bulk supersymmetry is similar to the SM case. The GUT scale 
values are derived using the supersymmetric RGE at the two loop instead of the SM ones. For 
the neutrinos however, one loop RGE were used with experimental inputs at the weak scale. The 
running of the masses are not dependent on the mixing angles for a low tan/3. Supersymmetry 
threshold corrections can play an important role while deriving the running masses. Running 
masses in the supersymmetric framework were obtained using the relevant matching conditions. 
As is well known, these effects are significant at large tan/3 and the corrections to the neutrino 
running through Yd and Y^ were considered [51]. 

The GUT scale masses and mixings chosen for the scan corresponded to tan/3 = 10 and are 
given in Table[VI] and [VII]. The results of the the scan i.e, the ranges for the c parameters are 
weakly dependent on tan/3 and can be applied for studying phenomenology for up to tan/3 ~ 25. 



TABLE VI: GUT scale Masses with supersymmetry for ta.n[3 = 10 



Mass 

(MeV ] 



Mass 

(GeV) 



Mass 
MeV 



Mass squared Differences 

eU2 



0.49+°;20 m. 



0.2361^:;^^^ TOe = 0.28 



+0.0000007 
-0.0000007 



TUd = 0.70 



+0.31 
-0.31 



?n, = 13 



+4 



TTib ^ 0.79 



-0.04 
-0.04 



59.9 



-0.000005 
-0.000005 



Aml^ = 1.6i0J° X 10-4 
Am|3 = 3.2+°;^3 x lO^^ 



nit = 92.21? j 



nir = 102ltJ5;^i 
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TABLE VII: Mixing angles for the quarks and leptons at GUT scale with supersynimetry for tan/3 = 10 



mixing angles (CKM) 



Mixing angles (PMNS) 



n 99«+0. 00087 



723 



= 0.0415: 



-0.00019 



13 = 0.0035t°;Sgl 



0.59 



+0.02 
-0.015 



^23 = 0.79tl]:i2 



'13 



0.154 



-1-0.016 
-0.016 



1. Quark Case 

The range chosen for the scan are the same as that for the SM case i.e. —2 < cq^^q^ < 4, 
—3 < CQg < 1 for the doublets. —2 < Cd^^d2,d-i < 3.5, for the down type singlets and —2 < Cui,u2 < 4, 
—4 < Cu3 < 1 for the up type singlets. The regions of c parameter space which satisfy the constraint 
of < x^ < 10 for the chosen scanning range are shown in Fig. (4) in Appendix B and the ranges 
are outlined in Table [VIII]. 

TABLE VIII: Ranges for the scanned regions of bulk hadronic parameters which satisfy < x^ < 10 for the 
supersymmetric case. 



parameter 


range 


parameter 


range 


parameter 


range 


CQs 


[-0.16,3.12] 
[-1.32,2.34] 

[-3,1] 


CD2 

CDs 


[-0.5,4] 

[-1.9,2.5] 

[-2,1.7] 


cui 

CU2 


[-1.6,4.0] 

[-2,2.4] 
[-4,1.0] 



2. Leptonic case 

Similar to the SM scenario two cases of neutrino mass generation are considered. The GUT 
scale input values for the x^ is given in Table [VI] and [VII]. 

(a)LLHH case 

The results of the scan of the LLHH case is very similar for both the SM case and the supersym- 
metric case. The expression for the neutrino mass matrix is given in Eq.(12). For the neutrino 
sector we allow the 0(1) Yukawa coupling to vary between -10 and 10 with a minimum of 0.08 
while that for the charged leptons are varied between -4 and 4 with a minimum of 0.08. The 
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doublets were scanned between -1.5 and 0.5 while the charged singlets were scanned between and 
4. The ranges for the c parameters for the LLHH case for the chosen scanning range satisfying the 
constraint < x^ < 10, is presented in Table[IX] and the plots for the c values are presented in 
Figs. (5) in Appendix B. 

TABLE IX: Ranges for scanned regions of the bulk leptonic parameters for the LLHH scenario in the 
supersymmetric case which satisfy < x^ < 10. 



parameter 


range 


parameter 


range 


CLi 


[-1.5,-0.22] 


CBi 


[2.6,3.7] 


CL2 


[-1.5,0.08] 


CE2 


[2.0,2.57] 


CL3 


[-1.5,0.04] 


C£3 


[1.1,1.8] 



(b)Dirac Neutrinos 

The expression for the mass matrix for the all the leptons is given by Eq.(lO) The scanning 
range for the c values of all the doublets and charged lepton singlets was in the range -1 to 
4.5, while the neutrino singlets were scanned in the range 3.5 to 9. The magnitude of 0(1) 
Yukawa parameters were varied between 0.08 and 4. The regions of the c parameters satisfying the 
constraint < x^ < 10 for the scanned ranges are presented in Table[X]. The ranges are presented 
in Fig. (6) in Appendix B. 

TABLE X: Ranges for the scanned regions of the bulk leptonic parameters for the Dirac case with super- 
symmetry which satisfy < x^ < 10 for the supersymmetric case. 



parameter 


range 


parameter 


range 


parameter 


range 


Cii 


[-1,2.6] 


CE^ 


[-0.86,3.46] 


CNi 


[5.68,8.9] 


CL2 


[-0.99,2.21] 


CE2 


[-1,2.24] 


CN2 


[5.67,8.99] 


CL3 


[-1,1.54] 


CEs 


[-1,1.49] 


CN3 


[5.64,8.99] 



To summarize, on comparing the SM and the SUSY fits, we find that within a given generation, 
the fields have a tendency to be localized slightly towards the IR for the SUSY case than for the 
SM case. This effect is more pronounced in the down sector and increases with tan/3. A comparison 
between the fits for the SM case and the SUSY case for tan/3 = 10 and 50 are presented in Figs. [7] 
in Appendix C. The underlying features of the fit in which the first two generations including the 
neutrinos are elementary from the ADS/CFT point of view while the third generation fermions 
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{tLiiR) having a tendency to be partially composite or composite, is maintained for both the SM 
and the SUSY case. From the choice of the c parameters, we find that the LLHH case admits a 
better fit to the neutrino data than the Dirac case. This is contrary to the observations made in 
[30] in normal RS where the c parameters for all the leptons were close to unity. It thus offered a 
more viable alternative than the LLHH case. In order to compensate for the weak warp factor in 
the Dirac case the right handed neutrinos had bulk masses cn ~ 7. This weak warping however, 
works in favour of the LLHH case where for c < 0.5 the effective 4D suppression scale is of the 
0{Mgut) resulting in fits with c parameters closer to unity. 

IV. SUSY SPECTRUM AND FLAVOUR PHENOMENOLOGY 

There are several ways to break supersymmetry within this RS set up at the GUT scale (see 
for example discussion in [31, 34-36, 52]. In the present work, we will consider only one particular 
set up which manifestly demonstrates the flavour structure of the fermions within the soft terms. 
More detailed analysis of supersymmetric spectrum will be addressed in [53]. We will assume in 
the following that supersymmetric breaking happens on the IR brane, or the GUT brane. Unlike 
the work of [33] and [36] we will not arrange the SM fields in any particular GUT representation. 
As has been discussed in these works, a GUT structure can be arranged with possible solutions for 
proton decay and doublet triplet splitting. Instead we parameterize SUSY breaking in terms of a 
single four dimensional spurion chiral superfield, X = 9'^F, which is localized on the GUT brane. 
However, for soft masses generated at the Planck brane as in [31], one may then impose the GS 
anomaly cancellation conditions on bulk masses to ensure unification of couplings at the Planck 
scale. We do not impose any such conditions as the soft masses are generated at the GUT scale. 

In the limit, the higher KK modes are decoupled from the GUT scale physics [38], the Kahler 
potential relevant for the scalar mass terms is given by 

/C(^) = / dy6{y-7rR)e-^''^'^k-^X^X [p,,^jQlQj + l3u,ijUJUj + |3d,^JDlDj + -fu,dHlaHu,d + ■■) 

(15) 

where /3 have dimensional carrying negative mass dimensions of -1 ( as the matter fields are five 
dimensional). ju,d s^-ie 0(1) parameters. 

The sfermion mass matrix is generated when the X fields get a vacuum expectation value 



m| r^ k ^ < X >' < X > Q'Q. The mass matrix will however not be diagonal in flavour space. 
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In the canonical basis, (15), the mass matrices take the form 

{^}h = ^1/2 k e(i-^-^^)^^-e(c.)e(c,) (16) 

where /3jj = 2kf3ij are dimensionless 0{1) parameters, ^(cj) are defined in Eq.(4). And the gravitino 
mass is defined as 



"^3/2 = -^^ = -Jar (17) 



The Higgs fields are localised on the GUT brane, their masses are given by m\j ^ = ^u4 "T'3/2- 

The A-terms are generated from the higher dimensional operators in the super potential of the 
type : 

VF(^) = j dy5{y - T:R)e-^^yk-^X [A^t^HuQ^Uj + Af^HaQidj + A^^HaLiEj + . . .) (18) 

where the A are dimensionful parameters having mass dimension -1. Substituting for the vev of 
the X, we have for the four dimensional trilinear couplings at the GUT scale: 

A^f = my,A[/'-^^-^'.^'''-ac^W,) (19) 

where we defined the dimensionless 0(1) parameters as A' = 2kA. The structure of the A terms 
and the corresponding fermion mass matrix are similar and they differ only by the choice of the 
0(1) parameters. Choosing A' = 2kY' , makes the down sector A terms diagonal in the mass basis 
of the fermions at the GUT scale. Henceforth, we shall work in this basis, with the 0(1) parameters 
of the A terms proportional to the 0(1) Yukawa parameters. 

The masses for the gauginos are obtained from the following operator in the lagrangian 

C= f d^ek'^XWAa'^A (20) 

At Mgut their masses will be be w-x /2 = fm^j2 where / is a 0(1) parameter, m^ /2 will be treated 
as an independent parameter. They are independent of the position of localization of X. as the 
profile for the gauginos is fiat corresponding to a bulk mass parameter of 0.5 [38, 40]. 

While the above equations set the boundary conditions at the high scale, the weak scale spec- 
trum is determined by the RGE evolution. In the present case, the spectrum at the high scale 
is completely non-universal as determined by the profiles of the zero modes of the matter chiral 
superfields. The structure of soft terms discussed here is similar to the ideas of flavourful super- 
symmetry discussed by [54, 55] and more recently by [56]. In the following we will present two 
example points one for the LHLH higher dimensional operator case and another for the Dirac case. 
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To begin with, in both the examples, we consider that all the 0(1) parameters appearing in the 
definitions of the soft parameters are proportional to the unit matrix. We will explicitly mention 
any deviations as required by the phenomenology when presenting numerical examples. This would 
mean that the matrices, A' , /3 in Eqs. (16, 19) are proportional to unit matrix and the parameters 
")u-,ld in Eq.(15) are equal to one. However, as we will see below, they play an important role 
in low energy phenomenology and one might frequently require to vary them within the 0(1) 
range, to satisfy phenomenological constraints. While studying the flavour phenomenology, we 
make sure that the soft terms are present in the super-CKM basis. The low-energy spectrum has 
been computed numerically using the spectrum generator SUSEFLAV [57]. 

3. LHLH operator case 

In this case we consider the following point, (21), in the c parameter space. It has a x^ of 5.5 
for the hadronic sector and 0.7341 for the leptonic sector. As expected it has a mostly composite 
right handed top quark. In addition, the leptonic doublets are also significantly composite in this 
case. 



cq^ = 2.740 CD, = 0.722 cu^ = 0.4024 cl, = -1.497 ce^ = 3.634 
CQ^ = 1.920 CD2 = 0.729 cu^ = 0.0652 cl^ = -0.224 ce^ = 2.290 
CQ3 = 0.960 CD3 = 0.801 ct/g = -3.5615 cl^ = -1.0738 ce,, = 1.769 (21) 

The choice of 0(1) parameters in the soft sector plays a role in determining the nature of the low 
energy spectrum. For a given set of c parameters, a naive choice of one for all the 0(1) parameters 
in the soft sector may or may not lead to an acceptable spectrum at Mgusy For the LLHH case 
corresponding to the choice in Eq.(21) the 0(1) parameters for all the soft masses are taken to 
be 1. The 0(1) parameters for the A terms are chosen to be A^ = 1.02y " while A" = Y and 
A^ = 0.6y ^. Corresponding to these choices of the 0(1) parameters and the c values in Eq.(21), 
the soft breaking terms at the GUT scale in GeV are given as: 



mq 



At, 



0.001 -0.03 -0.27 
-0.03 0.85 7.6 
-0.27 7.6 68.7 

-0.002 -0.09 -0.84 

-0.0002 1.09 -6.3 

-10-*^ 0.01 439.4 



\mu 



■.A 



D 



11.5 -63.5 156.1 

-63.5 349.6 -859.1 

156.1 -859.1 2110.8 

0.03 

-0.40 

-40.6 



;to_d = 



;mL 



105.03 -90.4 155.2 
-90.4 77.8 -133.6 
155.2 -133.6 229.5 

7.7 8.8 -147.5 

8.8 9.9 -167.0 
-147.5 -167.0 2798.4 
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ruE 



0.001 -0.017 0.05 

-0.017 0.26 -0.89 

0.05 -0.89 3.04 



]Ae = 



0.008 
1.81 
-31.3 



(22) 



A couple of interesting features of the above spectrum are (i) at least one of the soft masses 
is tachyonic (ii) significant amount of flavour violation present at the high scale. However at the 
weak scale, things are significantly different. This is because the RG running is quite different for 
the diagonal terms compared to the off-diagonal ones. In fact, the off-diagonal entries barely run, 
where as the corrections to the diagonal ones are quite significant. As an illustration, consider the 
slepton mass matrix at the weak scale. The analytic output at the weak scale for the diagonal 
terms can be approximated as 






(23) 



mi^{m 



susyj 



which receive gauge contributions while the off diagonal elements do not. The A terms are not 
large enough to make the off diagonal elements of the soft mass matrices comparable with the 
diagonal terms. An example for the sleptons for the case under consideration is given as 

2.2 X 10^ 77.0 -2.1 X 10^ 

77.0 2.2 X 10^ -2.7 x 10^ GeV^ 

-2.1 X 10"^ -2.7 X 10"^ 7.7 x 10^ 

1.1 X 10^ 3.4 X 10"^ 2.2 X 10^1 
susyj - 3.4 X 10""^ 1.1 X 10^ 5.9 X 10^ GeV^ 
2.2 X 10-1 5.9 X 10^ 5.4 x 10^ 

We see that the off-diagonal entry has barely enhanced where as the diagonal entries have been 
significantly modified. Constraints from fiavour violation would restrict the mass scales of m^/2 and 
Ml/2- The most stringent constraints are from the transitions between the first two generations ie 
from K^ —7- K^ and /i — )■ e-l-7. The expressions for the Kl — Kg mass difference and the branching 
fractions for /i — )■ e + 7 can be found in [58, 59]. We impose the fiavour constraints from all the 
existing data on the 5 parameters. 



m^^m. 



(24) 
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TABLE XI: Experimental upper bounds on the (5''°"'" obtained for rh„ = 2.1 TeV and rhi = 0.7 TeV 



(iJ) 


\s?,\ 


\sL\ 


\SEr\ 


l-^fi^l 


Kl\ 


I-^IlI 


\sRr\ 


Kr\ 


12 


0.053 


0.0002 


0.0003 


3.8x10-6 


0.0003 


3.8x10-6 


0.03 


0.03 


13 


0.34 


0.14 


0.06 


0.03 


0.06 


0.03 


0.26 


- 


23 


0.61 


0.16 


0.01 


0.04 


0.02 


0.04 


0.84 


- 



Bounds from the flavour violating processes are obtained using the mass insertion approximation 
[58] and the results of [60] defined in the Super-CKM basis. The flavour violating indices are defined 

as Sij{i 7^ j) = ^^ — —{i 7^ j) are evaluated in the basis in which the down sector is diagonal, 

with U being the rotation matrix which rotates the corresponding fermion mass matrix. The S 
are evaluated at the weak scale and we scale the bounds of [60] to the present mass scales. In 
Table[XII] we present the low energy spectrum corresponding to the sample point in Eq.(21). The 
low energy 6's are presented in Table XIII. 

TABLE XII: Soft spectrum for LLHH case: nisusy = 1-06 TeV, rrig = 2.64 TeV, ^i = 3.43TeV, tan/S = 25 



Parameter 


Mass(TeV) 


Parameter 


Mass (TeV) 


Parameter 


Mass (TeV) 


Parameter 


Mass(Tev) 


Parameter 


Mass (TeV) 


ii 


0.47 


6i 


1.01 


Tl 


0.726 


i>T 


2.78 


^1 


0.465 


i2 


1.05 


&2 


2.14 


T2 


2.79 


/>M 


0.483 


^2 


0.929 


cr 


2.24 


SR 


2.40 


fj-R 


0.478 


Ve 


0.469 


A^3 


3.38 


cl 


2.48 


SL 


2.48 


ML 


1.05 


- 


- 


A^4 


3.39 


ur 


2.24 


Ar, 


2.40 


SR 


0.476 


- 


- 


Cl 


0.895 


ul 


2.48 


dL 


2.48 


sl 


1.05 


- 


- 


C2 


3.43 


ruAQ 


3.23 


m± 


3.23 


ruh 


0.12186 


rriH 


3.06 


- 


- 



TABLE XIII: Low energy 5' s for quarks and leptons corresponding to the points in Eq.(21) for the LLHH 
case evaluated for fhq ^ 2.1TeV and m; = 0.7 TeV 



(iJ) 


l'^?.! 


I^^lI 


\SEr\ 


Kr\ 


\SRl\ 


Kl\ 


Kr\ 


Kr\ 


Kr\ 


12 


0.0003 


0.0001 


10-1" 


10-8 


10-8 


10-5 


0.001 


10-10 


0.001 


13 


0.01 


0.04 


io-« 


io-« 


10-5 


0.01 


0.005 


10-^ 


0.01 


23 


0.05 


0.05 


10-6 


10-5 


10-4 


0.04 


0.003 


0.0001 


0.07 
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4- Dirac Case 

For the case where neutrinos are of Dirac type the c parameters in Eq.(25) with x^ of 0.3211 for 
the hadronic sector and 0.1481 for the leptonic sector were chosen. The c values for the doublets 
in this case indicate they are predominantly elementary from the CFT point of view especially for 
the first two generations. The third generation however may be partially composite as in this case. 

CQ^ = 1.895 CD, = 1.898 cu, = 1.738 cl, = 1.293 ce, = 2.480 ctv, = 6.783 
CQ^ = 1.467 CD2 = 1-271 CU2 = 1-124 cl^ = 1-311 ce^ = 1-406 cn^ = 7-346 
CQ3 = -0.137 CDs = 1-394 cus = -0.356 c^g = 0.260 ce., = 0.237 cn, = 7.332 (25) 

The generic feature of soft mass matrices, discussed in the LLHH case, of the spectrum being 
tachyonic at the high scale and the diagonal terms evolving more than the off diagonal elements 
apply to this case as well. Corresponding to the c values in Eq.(25), the soft masses at the GUT 
scale have been evaluated for rn^V^ = 800 GeV while choosing Mi^ = 1200 GeV for the three 
gauginos. The 0(1) parameters corresponding to (w,q)33 and (m[/)33 were chosen to be 4 while for 
the others they are set to be 1. A^ = 1.151^ '^\/i,j, for all the A terms of the up sector while for the 
down sector and the leptons they were set equal to the corresponding 0(1) Yukawa couplings. The 
high scale soft breaking matrices in GeV are given in Eq.(26). In Table[XIV] we present the low 
energy spectrum corresponding to the sample point in Eq.(25). The low energy 5' s are presented 
in Table XV. 



mq 



0.60 1.2 30.0 
1.2 8.7 -19.4 
30.0 -19.4 2560.6 



■,mu 



0.7 -3.8 17.1 
-3.8 32.1 70.3 
17.1 70.3 2971.6 



■,mD 



0.10 -0.85 -1.7 
-0.85 7.1 14.6 
-1.7 14.6 29.8 



A 



u 



lO"-'^ -0.42 -5.8 
lO"-'^ 1.18 0.20 
10-5 -0.005 488.8 



:A 



D 



0.03 

0.57 
-40.6 



■,rnL 



8.7 0.9 62.5 
0.9 0.1 7.0 
62.5 7.0 446.3 



niE 



0.2 -0.4 10.5 

-0.4 0.7 -17.9 

1.0 -17.29 442.7 



■,Ae 



-0.01 

-3.01 
52.1 



(26) 
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TABLE XIV: Soft spectrum for Dirac case: rususy = 1-05 TeV, rrig = 2.65 TeV, ^ = 4.32TeV, tan[3 = 25 



Parameter 


Mass(TeV) 


Parameter 


Mass(TeV) 


Parameter 


Mass(TeV) 


Parameter 


Mass(Tev) 


Parameter 


Mass (TeV) 


ii 


0.702 


6i 


2.06 


Tl 


0.480 


Dr 


0.570 


^1 


0.465 


h 


2.31 


h 


2.32 


T2 


0.802 


/>M 


0.624 


N^ 


0.928 


cr 


2.25 


SR 


2.36 


fj-R 


0.608 


Ve 


0.625 


N3 


4.26 


cl 


2.45 


SL 


2.45 


Al 


0.902 


- 


- 


7V4 


4.26 


UR 


2.25 


Ar 


2.36 


SR 


0.610 


- 


- 


Cl 


0.894 


ul 


2.45 


dL 


2.45 


Cl 


0.903 


- 


- 


C2 


4.32 


niAQ 


4.18 


m± 


4.18 


mh 


0.1235 


niH 


3.96 


- 


- 



TABLE XV: Low energy S's for the Dirac Case corresponding to the point in Eq.(25) for Dirac case evaluated 
for rhq = 2.1TcV and m; = 0.7 TeV 



(ij) 


\s2,\ 


\sL\ 


\SSr\ 


\sU 


Kl\ 


Kl\ 


\5Rr\ 


1-51^1 


Kr\ 


12 


0.0003 


10^6 


10-10 


10-*^ 


10-s 


10-5 


10-7 


10-7 


0.00005 


13 


0.01 


0.007 


io-« 


10-s 


10-5 


0.002 


10-6 


10-4 


0.06 


23 


0.06 


10-4 


10-6 


10-5 


10-5 


0.01 


10-4 


0.0006 


0.001 



V. OUTLOOK 

The Randall-Sundrum framework is typically considered to be the geometric avatar of the 
Proggatt-Nielsen models. In the present work, we have considered a warped extra dimension close 
to the GUT scale. We fit the quark masses and the CKM mixing angles and determined the range 
of the c parameters which give a reasonable x^ fit. The 0(1) parameters associated with the 
Yukawa couplings have also been varied accordingly. Though the top quark Yukawa is smaller at 
the high scale compared to the weak scale, it is still large enough that one requires a large negative 
bulk mass parameter for the right handed top quark. For the leptons, we considered two particular 
models for neutrino masses (a) with Planck scale lepton number violating operator and (b) Dirac 
neutrino masses. 

The results show that there is a significant difference in the RS models at the weak scale and 
the RS models at the GUT scale especially if one focuses on the neutrino sector. In the weak scale 
models, the Planck scale lepton number violating higher dimensional operator was very hard to 
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accommodate with perturbative Yukawa couplings and thus was highly disfavoured. The Dirac and 
the Majorana cases were favoured though they were strongly constrained by the data from flavour 
violating rare decays. The situation is sort of reversed in the GUT scale RS models, though not 
exactly. The higher dimensional Planck scale operator fits the data very well, where as the Dirac 
case requires larger c values, some times with bulk mass parameters almost an order of magnitude 
larger than the cut-off scale. In the hadronic sector, the situation is not so dramatic. The top 
quark Yukawa though reduces at the GUT scale compared to the weak scale, still requires that the 
right handed top to be located close to the IR, making it a composite as in the weak scale models. 
The supersymmetric version of the same set up is far more interesting as it can have possible 
observable signatures at the weak scale. The main difference in fitting of the fermion masses is due 
to the presence of the additional parameter tan/? in the supersymmetric case. The difference in the 
c values is not very significant for low values of tan (3. At large tan /3, for a given generation, the zero 
modes are localized more towards the IR brane as compared to the SM case. This effect is more 
pronounced in the down sector as shown in Figures [7]. We parameterise SUSY breaking by a single 
spurion field localised on the IR brane. The resultant soft masses depend on the profiles of the zero 
modes of the chiral superfields and contain fiavour violation. At the weak scale the constraints 
from first two generation fiavour transitions rule out light spectrum. Another significant constraint 
comes the light Higgs mass at 125 GeV. The trilinear couplings have the same form as the Yukawa 
couplings in this model as long as the 0(1) parameters associated with both the parameters are 
taken to be proportional to each other. If all the 0(1) parameters at high scale are take to be 
exactly unity, the weak scale values of At are small to generate a 125 GeV light Higgs, for stops 
of masses ~ 1.5 — 2 TeV. However, with a minor variation of the 0(1) parameters for At at the 
high scale, 125 GeV Higgs is easily possible. More variations of supersymmetric breaking and the 
corresponding spectra will be discussed later [53]. 
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APPENDIX A: PLOTS FOR RANGES OF C PARAMETERS FOR QUARKS AND LEP- 
TONS FOR SM FITS AT THE GUT SCALE 

1. Range of c parameters for the quarks 



a a 
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C 







^ ■ 


■ . ■'•.■■'^". 

■ I . ' 

m • r ^ 





^mm 




■ 



'. '. :;-'A' ... 






'sW^ff^ 



ij 











1' 

: . . ■ '^] 




1 , , U , ■ 1 ■ ' 1 !■-! 


1 1 1 1 

■..;■■. .'■■ 
*?r ■ J** 1 







FIG. 1: The points in the above figures correspond to a x^ between 1 and 10 for the SM fits at the GUT 
scale. The plot represents the parameter space for the bulk masses of the quarks. 
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2. Range of c parameters for the for leptons for the LLHH scenario 









- 



m 











m -*-r.r".- !- . 



-""^r^^ _ 



^ 




= 1-2 





-w 


- 


- 



■IJ 

--E2 



FIG. 2: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the LLHH case. 
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3. Range of c parameters for the leptons. 
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FIG. 3: The points in the above figures correspond to a x^ between 1 and 10 for the SM fits 
scale. The plot represents the parameter space for the bulk masses of the leptons corresponding 
type neutrinos. 
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APPENDIX B: PLOTS FOR RANGES OF C PARAMETERS FOR QUARKS AND LEP- 
TONS FOR THE SUPERSYMMETRIC CASE 

1. Range of c parameters for the quarks 
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FIG. 4: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the quarks. 
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2. Range of c parameters for the for leptons for the LLHH scenario 
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FIG. 5: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the LLHH case. 
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3. Range of c parameters for the for leptons for the case of Dirac neutrinos. 
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FIG. 6: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the Dirac type neutrinos. 
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APPENDIX C: COMPARATIVE PLOT BETWEEN SM AND SUPERSYMMETRIC 
FITS FOR FERMIONS 
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FIG. 7: The plots correspond to eomparison between SM and supersymmetric fits for down sector fermions. 
Red corresponds to SM, wliile Blue and Green corresponds to supersymmetric case for tan/3 = 10 and 
tan/3 — 50 respectively 
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Warped Alternatives to Proggatt-Nielsen Models 
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Abstract 
We consider the Randall-Sundrum (RS) set-up to be a theory of flavour, as an alternative to Froggatt- 
Nielsen (FN) models instead of as a solution to the hierarchy problem. The RS framework is modified by 
taking the low energy brane to be at the GUT scale. This also alleviates constraints from flavour physics. 
^Z Fermion masses and mixing angles are fit at the GUT scale. The ranges of the bulk mass parameters are 

^H determined using a y^ fit taking in to consideration the variation in 0(1) parameters. In the hadronic 

O , sector, the heavy top quark requires large bulk mass parameters localising the right handed top quark close 

'^N to the IR brane. Two cases of neutrino masses are considered (a) Planck scale lepton number violation and 

LJ (b) Dirac neutrino masses. Contrary to the case of weak scale RS models, both these cases give reasonable 

^_^ fits to the data, with the Planck scale lepton number violation fitting slightly better compared to the Dirac 

case. In the Supersymmetric version, the fits are not significantly different except for the variation in tan/3. 
If the Higgs superfields and the SUSY breaking spurion are localized on the same brane then the structure 
of the sfermion masses are determined by the profiles of the zero modes of the hypcrmultiplcts in the bulk. 
^^ Trilinear terms have the same structure as the Yukawa matrices. The resultant squark spectrum is around 

^^^J ^2 — 3 TeV required by the light Higgs mass to be around 125 GeV and to satisfy the flavour violating 

OO constraints. 
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I. INTRODUCTION 

One of the celebrated solutions of the fermion flavour problem is the Proggatt-Nielsen Mechanism 
[1]. According to this prescription, the symmetry group of the Standard Model (SM) is augmented 
by a horizontal U{l)x group under which all the SM fermions and the Higgs field are charged. The 
effective theory includes a flavon field X and the Yukawa couplings are generated from the higher 
dimensional operators which are invariant under the U{l)x and the Standard Model (SM) gauge 
group. For example, the up-type quark mass matrix has the form: Yj^lj^j^^^^ '^"^ '^""QiH^Uj, 
where c/ is U{l)x charge of the / field and i,j are the generation indices. The flavon field X 
develops a vacuum expectation value (vev) such that 0.22 ~ Ac ~< X > /Mpi, Ac being the 
Cabibbo angle. Y-^ are taken to be 0(1) parameters. Fermion mass matrices including their 
mixing patterns can be fit to the data by choosing appropriate U{l)x charges for various fields. 
An UV completion of the model can be constructed by including heavy chiral fermions in to the 
theory; integrating these heavy fields would lead to the relevant non-renormalizable operators (for 
a review, see [2] ). 

The U{l)x symmetry introduces additional anomalies in to the theory and subsequently, strong 
constraints on the U{l)x charges for various fields. In supersymmetric models with a single flavon 
fleld, one typically has to resort to Green-Schwarz (GS) mechanism to cancel the anomalies. The 
solution set of U{l)x charges for the fermions and the Higgs which satisfy the fermion data as well 
as the anomaly cancellation^ requirement have been studied in [3-14] and recently updated in [15]. 
These models typically lead to large flavour violations at the weak scale in gravity mediated super- 
symmetry breaking models due to contributions from the U{l)x D-terms. While the constraints 
from the flavour sector on the available solutions are very tight, it may still be possible to ease 
them without requiring the superpartner masses to be very high [16, 17]. The flavour constraints 
may also be alleviated to some extent by considering U{1) x U'{1) class of models [18, 19]. In 
the present work, we will study the extra- dimensional alternative [20] to understand the flavour 
hierarchy in particular concentrating on the supersymmteric Randall- Sundrum (RS) set up. 

The Randall-Sundrum framework [21] which elegantly provides a solution to the hierarchy 
problem via warping in the extra dimensional space can also thought to be a theory of flavour. 
It has been observed sometime ago that the flavour changing neutral currents (FCNC) can be 
suppressed due to the so-called RS-GIM mechanism [22]. However, in the absence of additional 



^ However, with two singlet fiavons there exist a unique solution which is completely non anomalous [3] 



flavour symmetries the constraints from FCNC are still very strong( [23-27]) (Detailed analysis for 
the hadronic sector can be found in [28, 29] and references there in. For a recent thorough analysis 
in the leptonic sector, please see [30]). Given these strong constraints on the RS set up at the weak 
scale, one can ask the question whether RS is suitable to be a theory of flavour as well as a solution 
to the hierarchy problem simultaneously. It might be that RS as a theory of a flavour might be 
better suited at the GUT scale rather than at the weak scale. The Froggatt-Nielsen models are 
typically defined at scales closer to the Planck scale, so perhaps flavour physics might have its 
origins at the Planck scale. 

With this rationale, in the present work we will consider RS to span between the Planck and 
the GUT scales. The fermion masses are flt in terms of the bulk mass parameters of the various 
fields, which take the role of the U{l)x charges of the FN mechanism. However, these parameters 
are less constrained compared to the U{l)x charges, as no additional conditions such as anomaly 
cancellations are required on them. While this has been the common understanding, in [31] it was 
pointed out that imposing unification conditions on gauge couplings in a theory with localization 
of fermions or hierarchical wave functions leads to strong constraints which are exactly in the same 
way as the Green-Schwarz anomaly cancellation conditions [32] ^ In the present setup we do not 
impose these conditions. 

Extra dimensions at GUT scale were considered in [33] while the RS version was considered 
by the authors in [34, 35] and later by the authors in [31]. Our work, however, is very closely 
related to the work of [36] who have done a thorough analysis of fermion mass spectrum, weak 
scale supersymmetric spectrum and fiavour phenomenology, assuming a particular Grand Unified 
Theory (GUT) model in such a RS setting. However, differences exist. In the present work we 
have not assumed any specific GUT model. Furthermore, we have used a frequentist approach 
to do the fermion mass fitting. While this makes it hard to directly compare the results between 
the two works, we hope they provide a complementary set of results. We also have taken in to 
consideration the constraints from neutrino masses and mixing angles which can have a significant 
effect on the lepton fiavour violation and slepton decays. 

The equivalent description for the RS set up in four dimensions can be thought of as a composite 
Higgs coupled to fermions with couplings which parameterise the 'partial compositeness' of the 
fermions[37]. In SUSY case, this partial compositeness can also affect the structure of the soft 



^ Typically applied in FN models, the Green-Schwarz anomaly cancellation conditions requires the anomaly factors 
to be in a particular ratio such that they are cancelled in String theory. 



masses. 

In the first part of our work, our aim lias been to provide a range of bulk mass parameters which 
fit the fermion masses and mixing patterns at the GUT scale. We believe this can be useful for 
model builders and other phenomenologists working in flavour physics and looking for an alternative 
to FN models. We have considered both supersymmetric as well as non supersymmetric versions 
of the RS framework at the GUT scale while fitting the data. The supersymmetric case has the 
added advantage that it could lead to observable signatures at the weak scale. We consider the 
case where SUSY breaking is considered to be on the same brane as where the Higgs is localized, 
which is the GUT brane. In this case, the sfermion mass matrices are determined by the zero mode 
profiles of the corresponding N=l superfields and thus the information of the fermion masses is 
propagated in to the soft sector. It is far more striking for the A-terms which follow the same 
structure as the Yukawa couplings. The spectrum is highly non-universal at the high scale, but, 
its pattern is constrained due to the ranges of bulk mass parameters which are in turn are fixed by 
their fits to fermion masses. The running effects make the diagonal terms large at the weak scale. 

The rest of the paper is organized as follows. In section II, we detail the RS setup we consider 
and derive the structure of the fermion masses. In section III we present the fermion mass fits 
and present the ranges for the bulk mass parameters for both the non-supersymmetric and the 
supersymmetric cases. In section IV we address the issue of supersymmetric breaking and derive 
supersymmetric spectrum for a particular supersymmetric breaking case. We end with summary 
and outlook in the last section. In Appendices A , Band C, we have presented plots relevant for 
fermion mass fits. 

II. RS AS A THEORY OF FLAVOUR 

The Randall-Sundrum frame work consists of two branes separated by an single warped extra 
dimension [21]. The line element for the RS background is given as 

ds^ = e-'^'^yri^^dx^dx'' - dy^ (1) 

where < y < ttR. Here y = is identified as the position of the UV brane and y = irR is the 
position of the IR brane. The scale associated with physics on the UV brane is Mpi while that on 
the IR brane is TeV. The solution to the hierarchy problem is achieved by exponential warping of 
scales i.e, Mpi = e^^"" M^cak, where kR ~ 0(11). 

In the modified set up we consider here, the scale associated with the IR brane is Mqut- This 



can be achieved by choosing kR ~ 1.5. We define the hierarchy between the scales, e, for this 
scenario to be 



10~^ (2) 



Mpianck 

We consider both supersymmetric and non-supersymmetric matter fields to propagate in the bulk. 
Localisation of the respective zero modes is dependent on the corresponding bulk masses. In both 
the cases, we assume that the Higgs field (two Higgs fields in the case of supersymmetric models) are 
localised on the GUT brane. We now proceed to briefly review the derivation of the mass matrices 
and their dependence on the zero mode profiles in both supersymmetric and non-supersymmetric 
cases. A couple of points are important to note at this juncture. Firstly, the typical lowest KK 
mass for a warped background is given as rriKK = e~^^'^k. In the present set up, the 'large' 
warp factor ensures the lowest KK modes are very heavy i.e, rnxK = eA^ ~ Mqut and thus are 
decoupled from low energy phenomenology. We do not consider their effects in this work for low 
energy phenomenology. Secondly, it turns out that the dependence of the zero mode mass matrices 
on the profiles is very similar in supersymmetric and non-supersymmetric cases. However, the 
fermion mass data at the high scale in the supersymmetric case could be different from the SM 
one, due to the dependence on tan/3 as well as the different RGE for the Yukawa coupling as we 
will discuss in the next section. 

A. Standard Model case 

The non-supersymmetric case or the Standard Model case has been studied in many works 
[23, 28-30]. The main difference in the present case is that while those studies have considered a 
kR ~ 0(11), while in the present case it is 0(1). We thus assume a grand desert from the weak 
scale to the GUT scale, where RS framework sets in. No attempt is made to solve the hierarchy 
problem, but the fiavour problem has a solution in terms of the localisation of the fields in an extra 
dimension at the GUT scale. We follow the notation of [30] and present the final formulae for 
the Yukawa mass matrices. The details of the KK expansion and the corresponding ortho-normal 



relations can be found in [301 and references therein. The five dimensional action has the form: 



'-' — '-'kill + OYuk -\- Ju ~^ ^higgs 

Skin = j d'^x j dy ^/^ {L{ilJ)-mL)L + E{ilt)-mE)E + ... ) 
5Yuk = I d'^x I dy 7=^ ( YuQUH + YdQDH + YeLEH^ 5{y - ttR) 
Sy = f d'^x f dy ^/^(-^LHLH or YnLNh) 5{y - ttR) (3) 

where we used the standard notation with the Q, U, D standing for the quark doublets, up-type and 
down type singles respectively, L and E, N stand for leptonic doublets and charged and neutral 
singlets respectively. H stands for the Higgs doublet with H = iaH*. We have suppressed the 
Higgs action in the above. Two specific ways for generating non-zero neutrino masses are considered 
(a) by a higher dimensional term localised at the GUT brane and (b) Dirac neutrino mass terms 
similar to the other fermions. A'^^ is the five dimensional reduced Planck scale ~ 2 x 10^® GeV. 

After the Kaluza-Klein (KK) reduction and imposing the orthonormal conditions, we can derive 
the 4D mass matrices for the zero-modes of the fermion fields. They have the form: 

where F stands for all the Yukawa matrices F = U,D,E and N, if the neutrinos have Dirac masses. 
Cj and c'- represent the bulk masses of the respective matter fields (second line of eq.(3)); defined 
as, for example, m^;- = C£;.A:. i,j denote the generation indices. If the neutrinos have Dirac masses 
then their mass matrix is given by Eq.(4). In case the neutrinos attain their masses through higher 
dimensional operator, the mass matrix is given by 

In Eqs. (4, 5), we have defined Y' = kY and k' = 2kK. These are dimensionless 0(1) parameters 
of anarchical nature^. Eq.(4) are used to fit all the fermion mass data at the GUT scale i.e, 
up and down type quark masses and the (Cabibbo-Kobayashi-Masakawa) CKM matrix, charged 
lepton masses, neutrino mass differences and the corresponding PMNS mixing matrix. In the case 
neutrinos get their masses through higher dimensional operator, Eq. (5) is used instead to fit their 
mass differences and mixing angles. 



Note that the Yukawa couphngs in Eq.(9) are dimensionful, with mass dimensions -1. 



B. Supersymmetric case 

In the supersymmetric case, the matter fermions are represented by hyper-multiplets^ propa- 
gating in the bulk. In terms of the 4D, N=l SUSY language, they can be expressed as two N=l 
chiral multiplets, $, ^'^. Following [38, 40], we write the 5D action in terms of two chiral fields with 
a (supersymmetric) bulk mass term to be 



85= f d^x I (fee-'^^y ($f <5 + ^^^"■f) + / Soe-'^^y^^ (dy + m$ - h^ $ 



(6) 



where M$ = c^k is the bulk mass. In writing the above, the radion field is suppressed by taking 
its vacuum expectation value, < Re[T) >= R. The super field <I>^ is taken to be odd under Z2. 
Thus, only <I> has a zero mode. Since we have a theory at the GUT scale, the KK modes can be 
considered to be decoupled from theory. In the effective theory, the profile of the zero mode of the 
<I> is determined by [38] 

5. - (^ - c) k^ /(°) = (7) 

Thus f^^' = e'a"'^' y . The superscript '"^ stands for the zero mode, which we will drop subse- 
quently^. In this effective theory, where the higher KK modes are completely decoupled, we can 
write the effective 4-D Kahler terms for the Z2 even zero modes as [31, 34, 35] 

where we have substituted for the profile solutions of (7). After integrating over the extra dimension 
y, the terms in Eq.(8) pick up a factor Zp = tt^ — \u {^^^~^ — l) where F = Q,U, D, L,E, 



(l-2cF)fe 



as 



before. We choose to work in a basis in which the Kahler terms are canonically normalized. We 
thus re-define the fields as $ — t- —r^=^. 

The effective four dimensional MSSM Yukawa couplings are determined from the superpotential 
terms written on the boundary. For Higgs localized on the IR brane, keeping only the zero modes 
of the chiral superfields, the effective four dimensional superpotential is given as [31, 38] 

+ e^l-^^^^'^ye^l-'^^^'yv^HoLiEj + . . ) 6iy - ttR) (9) 



'^ N=l Supersymmetry in 5D has the particle content of N=2 Supersymmetry in 4D. The hypers can be expressed in 

N=l, 4D language as two chiral superfields, where as the Vectors can be expressed as a vector and chiral superfield 

[38, 39]. 
"" In the component form, the scalar component and the fermion components of the chiral super field $ have different 

bulk masses. However, the solution for the profile for the scalar and the fermion components turns out to be the 

same. 



The Higgs fields are canonically normalized as H^^d ~^ ^^^^H^^d- In the canonical basis, after the 
fields have been redefined the fermion mass matrices can be derived from Eq. (9) to be 

{Mf\, = !^y/.e(i-^-^i)'=-C(c.)e(4) (10) 

where Ci,c', denote the bulk mass parameters for various fields, ^(cj) are defined in Eq.(4) The 
mass matrix in Eq.(lO) can be approximated as {MF)ij '^ ^7^C(l)e ~'^'^'^J' '"'^ where the ^(c) is 
absorbed into the 0(1) parameters Y' and is now collectively referred to as 0(1). This is true only 
as long as the c parameter lies between and 1. But as we have seen earlier, in some realistic cases 
especially related to neutrino masses and the top quark, the values of \c\ could be large to fit the 
data. Redefining the 0(1) Yukawa by absorbing the c parameters would shift the ranges of the c 
parameters far away from what they are, especially in the case where \c\ > 1. As in the SM case, 
we define dimensionless 0(1) Yukawa couplings as Y' = 2kY and are of anarchical nature. While 
the Dirac masses for the neutrinos have the same structure as the other fermion mass matrices, 
the higher dimensional operator has a different form determined by the super-potential term 



w(4) 



I'dydiy - 7ri?)e-3-fe) (^e^^l-^^^^^ye^l-'^J^^y ^HuHuLiLj') 



(11) 



The neutrino mass matrix in this case is given as 



where k' = 2kK is the dimensionless 0(1) parameters. The fermion mass matrices carry the 
same form as in the SM and the supersymmetric cases and thus their dependence on q and 0(1) 
parameters is the same. 



III. FERMION MASS FITS 

From the previous section, we have seen that in addition to the bulk mass parameters, the 0(1) 
Yukawa parameters also play a role in fixing the fermion masses and mixing angles. We fit the 
masses and the mixing angles of the quark sector and the neutrino sector at the GUT scale for 
both the SM and the supersymmetric cases. We will use a frequentist approach, i.e, we minimise 
the x^ function, which is defined as follows: 

^^ii f'-f' ^j (13) 



eory ^^ ^j^^ theory number for the j^^ observable and y^^^ 
quoted by experiments with a measurement uncertainty of aj. In the present case the theory 



where, y- ^"^^ is the theory number for the j^^ observable and y^^^ is its corresponding number 



parameters are just the bulk mass parameters and the 0(1) Yukawa entries in supersymmetric 
and non-supersymmetric cases. We define < x^ < 10 to be a good fit and we try to find regions 
in the parameters space of bulk mass parameters and 0(1) Yukawa parameters which satisfy this 
condition^. The 0(1) Yukawa parameters are varied between -4 and 4, with a lower bound of 0.08 
on \Y\ to avoid unnaturally small Yukawa parameters. As far as the bulk mass parameters are 
concerned, since they are given as ck, we prefer to vary the c parameters between —1 to 1, so as 
not to go beyond the 5D cut-off, k. This will remove any possible inconsistencies in the theory 
due to non-perturbative Yukawa couplings. However, as we will see it is not always possible to fit 
the data within this range of c parameters. We will mention the range chosen specifically for each 
case. 

The minimization of the x^ function was performed using MINUIT [41]. We can minimise 
the hadronic and the leptonic sectors independently as they are dependent on different sets of 
parameters which are uncorrelated. The methodology is similar to the ones used in fermion mass 
fitting in GUT models [42, 43] and also the one used in [30]. 

A. Standard Model (SM) Case 

In this section, we present the fits in the SM case. For the GUT scale values of the quark and 
lepton masses and CKM mixing matrices, we use the results of [44]. In the analysis of [44], two 
loop RGE have been used to run the Yukawa couplings of the up- type quarks, down- type quarks 
and charged leptons from the weak scale all the way up to the GUT scale. For the neutrino data 
we used the publicly available package REAP [45] to compute the high scale values. The masses of 
the SM fermions at the GUT scale used in our fits are presented in Table I. The CKM and PMNS 
mixing matrices are presented in Table II. 

1. SM Quark Sector fits 

The up and down mass matrices are given in terms of fermion mass matrix of Eq.(4). The 
theory parameters which are varied simultaneously to minimise the x^ ^^ Eq.(13) include: three 
CQ^, each of c„^ and c^^ and 18 0(1) Yukawa parameters. We would expect that the light quarks 



We will mention the results with lower x^ a-t relevant places. 



TABLE I: GUT scale masses of fcrmions for the SM case 



Mass 

(MeV ) 


Mass 
(GeV) 


Mass 
MeV 


Mass squared Differences 


TO„ = 0.48+0;?° 

ma = 1.14tl!i^ 

m, - 22+_l 


TO, = 0.235t°:°^^ 

TO, = i.ot°:Sl 
TOt = 74.01;^:? 


TO - n 4606+"°"°"°°°'' 
irie — u.^ioyo^Q gQQQooo4 

TO - qq i4+"°"°°°* 
in^ — ^^-^^-0.0000089 

rur = leSS.SStoig 


Amj^ = 1.5+°;2; X 10-4 
Am|3 = 4.6+°;^!^ X 10-3 



TABLE II: Mixing angles for the hadronic and the leptonic sector for the SM case 



mixing angles (CKM) 



^12 - U./;^b_o 00087 

Q — n nzLi c;+o.oooi9 
013 = 0.0035t°;l]°} 



Mixing angles (PMNS) 



0.59: 



0.015 



023 = 0.79t^;i2 



'13 



= 0.154 



+0.016 
-0.016 



would be localised close to the UV brane ( c > 1/2 ) and the heavy quarks close to the IR brane 
(c < 1/2). However, for this particular range of 0(1) Yukawa parameters, it is difficult to fit the 
data for |c| within unity. We thus enlarged the range for the c parameters. 

The range chosen for the scan of the c parameters chosen is: —2 < cq-^^q.^ < 4, —3 < cqj < 1 for 
the doublets. —2 < Cd^^d243 < 3-^' fo'^ ^^^ down type singlets and —2 < Cu^^u2 < 4, —4 < 0^3 < 1 
for the up type singlets. We fit the quark masses and the CKM mixing angles at the GUT scale. 
The top quark is definitely lighter at the GUT scale, but still we see that most of the points that 
fit the data lie outside of |c| < 1. This is evident from the the negative values of the cqj and Cf/3 
that fit the data. 

The regions of c parameter space which satisfy the constraint of < x^ < 10 for the chosen 
scanning range are shown in Fig.(l) in Appendix A and the ranges are outlined in Table [III]. We 
see that the first two generation bulk mass parameters are concentrated on the positive c values 
where as the third generation, the doublet and more so the right handed top is localised close to 
the GUT scale brane. Comparing these results with that of the normal RS, we find that the masses 
for the light quark fields can be fit with c ^ 0.6 — 0.7. This can be attributed to the large warping 
where 0.5 < c < 1 is sufficient to reproduce the masses for light quarks [46]. 
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TABLE III: Allowed range of c parameters in the SM case. These parameters satisfy < x^ < 10 for the 
SM case. The corresponding figure is 1 in Appendix A. 



parameter 


range 


parameter 


range 


parameter 


range 


CQ3 


[0,3.0] 
[-1.95,2.36] 

[-3,1] 


CD2 


[0.78,4] 
[0.39,3.02] 
[0.39,2.21] 


CC/2 
CC/3 


[-0.97,3.98] 

[-1.99,2.43] 

[-4,1.0] 



2. SM Leptonic Mass fits 



Unlike the quark case, the fits in the leptonic sector are far more difficult and more constraining 
due to the small mass differences and the large mixing in the neutrino sector. As mentioned, we 
will consider two different cases of neutrino masses while fitting the leptonic data, 
(a) LLHH higher dimensional operator 

Planck scale lepton number violation is an interesting idea which manifests itself with higher 
dimensional operator suppressed by the Planck scale. In four dimensions such an operator generates 
too small neutrino mases. It is typically used as a perturbation over an existing neutrino mass 
model [47]. If not, it needs an enhancement of 0(10'^ — 10^) to be consistent with the data. In 
the standard RS framework close to the weak scale with bulk fermions, this higher dimensional 
operator is still constrained however for different reasons. While the neutrino masses can be fit by 
placing the doublet fields L close to the UV brane, the charged lepton masses become very tiny 
unless the singlet fields (E) are placed deep in the IR[30]. This leads to inconsistencies in the theory 
with large non-perturbative Yukawa couplings. The question arises whether the situation repeats 
itself when we consider the modified RS setup. This can be checked as follows. The neutrino 
masses are generated by the higher dimensional operator as given in Eq.(3). The corresponding 
neutrino mass matrix is given by Eq.(5) while the mass matrix for the charged leptons is given by 
Eq.(4). 



For simplicity assume cl- = c^V i. For cl < 0.5 the mass matrix in Eq.(5) becomes 



nil, = K 



2eA 



:i - 2cl) 



(14) 



It is clear that cl ~ —4 is required to get neutrino masses 0(0.04) eV for a warp factor for e ~ 10~^. 
As Cl increases, beyond 0.5, this formula is no longer valid, the neutrino masses become smaller 
and hence do not fit the neutrino mass data with 0(1) Yukawa couplings. Thus a mildly negative 
Cl should be able to fit the data without large inconsistencies. 
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A second enhancement can also come from the k', which is the corresponding 0(1) Yukawa. 
With this in mind, we enhance the range of the scanning of the Yukawa couphngs from 0.08 to 4 
to 0.08 to 10. This would help us to accommodate cl values close to ~ —1. The final scanning 
ranges we have chosen are: the doublets (c^-) are varied between -1.5 and 0.5, while the charged 
singlets were scanned between and 4. The region of c values which give a good fit to leptonic 
masses, i.e, satisfying the constraint < x^ < 10, is presented in Table[IV]. The plots for these 
ranges of c values are presented in Figs. (2) in Appendix A. 

TABLE IV: Ranges for scanned regions of the bulk leptonic parameters for the LLHH in the SM case which 
satisfy < x^ < 10. 



parameter 


range 


parameter 


range 


Cii 


[-1.5,-1.15] 


ce^ 


[2.8,4.0] 


CL2 


[-1.5,-0.97] 


CE2 


[1.8,2.4] 


CL3 


[-1.5,-1.22] 


CEs 


[1.2,1.69] 



(b)Dirac type Neutrinos 

The case of Dirac neutrinos is interesting possibility though it requires imposition of a global lepton 
number conservation^. The running of the neutrino masses from the weak scale to high scale is 
different in this case. However with a normal hierarchy of neutrinos and low tan/3 the differences 
are insignificant [49, 50]. 

Assuming that there is not much of a difference for normal hierarchy, we choose the following 
scanning range for the c parameters. The doublets {cl-) and charged lepton singlets (c^;-) are 
scanned within the range -1 to 4.5, while the neutrino singlets were scanned in the range 3.5 
to 9. Such a large value of the bulk mass parameters for the singlets is needed to suppress the 
corresponding neutrino masses sufficiently. The 0(1) Yukawa parameters were varied between 
0.08 and 4. Comparing the results of Dirac neutrino mass fits with that of the weak scale RS 
models, [30], we find that the cn are roughly a factor 7 — 8 larger compared to the cjsr at the weak 
scale. This is purely because of the weaker warp factor we are considering in the present case. 
Increasing the range of the 0(1) Yukawa parameters would only make things worse. The ranges 
for the c values corresponding to SM fits with Dirac neutrinos case are presented Table [V]. The 



''' In fact, it is possible to hide lepton number violation in this case through a careful location of the right handed 
fermion fields [48]. We will not consider this case here. 
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plots for the c parameters are presented case in Fig[3] in Appendix A. 



TABLE V: Ranges for the scanned regions of the bulk leptonic parameters for the Dirac case which satisfy 
< x^ < 10 for the SM case. 



parameter 


range 


parameter 


range 


parameter 


range 


CLi 


[ -1,2.9] 


CE^ 


[0.39,3.62] 


CATi 


[5.29,8.97] 


CL2 


[-0.99,2.7] 


CE2 


[-1.0,2.63] 


CN2 


[5.31,8.99] 


CL3 


[-0.99,1.98] 


CEs 


[-0.99,1.93] 


CNs 


[5.12,8.97] 



B. Supersymmetric Case 

The analysis for the case with bulk supersymmetry is similar to the SM case. The GUT scale 
values are derived using the supersymmetric RGE at the two loop instead of the SM ones. For 
the neutrinos however, one loop RGE were used with experimental inputs at the weak scale. The 
running of the masses are not dependent on the mixing angles for a low tan/3. Supersymmetry 
threshold corrections can play an important role while deriving the running masses. Running 
masses in the supersymmetric framework were obtained using the relevant matching conditions. 
As is well known, these effects are significant at large tan/3 and the corrections to the neutrino 
running through Yd and Y^ were considered [51]. 

The GUT scale masses and mixings chosen for the scan corresponded to tan/3 = 10 and are 
given in Table[VI] and [VII]. The results of the the scan i.e, the ranges for the c parameters are 
weakly dependent on tan/3 and can be applied for studying phenomenology for up to tan/3 ~ 25. 



TABLE VI: GUT scale Masses with supersymmetry for ta.n[3 = 10 



Mass 

(MeV ] 



Mass 

(GeV) 



Mass 
MeV 



Mass squared Differences 

eU2 



0.49+°;20 m. 



0.2361^:;^^^ TOe = 0.28 



+0.0000007 
-0.0000007 



TUd = 0.70 



+0.31 
-0.31 



?n, = 13 



+4 



TTib ^ 0.79 



-0.04 
-0.04 



59.9 



-0.000005 
-0.000005 



Aml^ = 1.6i0J° X 10-4 
Am|3 = 3.2+°;^3 x lO^^ 



nit = 92.21? j 



nir = 102ltJ5;^i 
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TABLE VII: Mixing angles for the quarks and leptons at GUT scale with supersynimetry for tan/3 = 10 



mixing angles (CKM) 



Mixing angles (PMNS) 



n 99«+0. 00087 



723 



= 0.0415: 



-0.00019 



13 = 0.0035t°;Sgl 



0.59 



+0.02 
-0.015 



^23 = 0.79tl]:i2 



'13 



0.154 



-1-0.016 
-0.016 



1. Quark Case 

The range chosen for the scan are the same as that for the SM case i.e. —2 < cq^^q^ < 4, 
—3 < CQg < 1 for the doublets. —2 < Cd^^d2,d-i < 3.5, for the down type singlets and —2 < Cui,u2 < 4, 
—4 < Cu3 < 1 for the up type singlets. The regions of c parameter space which satisfy the constraint 
of < x^ < 10 for the chosen scanning range are shown in Fig. (4) in Appendix B and the ranges 
are outlined in Table [VIII]. 

TABLE VIII: Ranges for the scanned regions of bulk hadronic parameters which satisfy < x^ < 10 for the 
supersymmetric case. 



parameter 


range 


parameter 


range 


parameter 


range 


CQs 


[-0.16,3.12] 
[-1.32,2.34] 

[-3,1] 


CD2 

CDs 


[-0.5,4] 

[-1.9,2.5] 

[-2,1.7] 


cui 

CU2 


[-1.6,4.0] 

[-2,2.4] 
[-4,1.0] 



2. Leptonic case 

Similar to the SM scenario two cases of neutrino mass generation are considered. The GUT 
scale input values for the x^ is given in Table [VI] and [VII]. 

(a)LLHH case 

The results of the scan of the LLHH case is very similar for both the SM case and the supersym- 
metric case. The expression for the neutrino mass matrix is given in Eq.(12). For the neutrino 
sector we allow the 0(1) Yukawa coupling to vary between -10 and 10 with a minimum of 0.08 
while that for the charged leptons are varied between -4 and 4 with a minimum of 0.08. The 
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doublets were scanned between -1.5 and 0.5 while the charged singlets were scanned between and 
4. The ranges for the c parameters for the LLHH case for the chosen scanning range satisfying the 
constraint < x^ < 10, is presented in Table[IX] and the plots for the c values are presented in 
Figs. (5) in Appendix B. 

TABLE IX: Ranges for scanned regions of the bulk leptonic parameters for the LLHH scenario in the 
supersymmetric case which satisfy < x^ < 10. 



parameter 


range 


parameter 


range 


CLi 


[-1.5,-0.22] 


CBi 


[2.6,3.7] 


CL2 


[-1.5,0.08] 


CE2 


[2.0,2.57] 


CL3 


[-1.5,0.04] 


C£3 


[1.1,1.8] 



(b)Dirac Neutrinos 

The expression for the mass matrix for the all the leptons is given by Eq.(lO) The scanning 
range for the c values of all the doublets and charged lepton singlets was in the range -1 to 
4.5, while the neutrino singlets were scanned in the range 3.5 to 9. The magnitude of 0(1) 
Yukawa parameters were varied between 0.08 and 4. The regions of the c parameters satisfying the 
constraint < x^ < 10 for the scanned ranges are presented in Table[X]. The ranges are presented 
in Fig. (6) in Appendix B. 

TABLE X: Ranges for the scanned regions of the bulk leptonic parameters for the Dirac case with super- 
symmetry which satisfy < x^ < 10 for the supersymmetric case. 



parameter 


range 


parameter 


range 


parameter 


range 


Cii 


[-1,2.6] 


CE^ 


[-0.86,3.46] 


CNi 


[5.68,8.9] 


CL2 


[-0.99,2.21] 


CE2 


[-1,2.24] 


CN2 


[5.67,8.99] 


CL3 


[-1,1.54] 


CEs 


[-1,1.49] 


CN3 


[5.64,8.99] 



To summarize, on comparing the SM and the SUSY fits, we find that within a given generation, 
the fields have a tendency to be localized slightly towards the IR for the SUSY case than for the 
SM case. This effect is more pronounced in the down sector and increases with tan/3. A comparison 
between the fits for the SM case and the SUSY case for tan/3 = 10 and 50 are presented in Figs. [7] 
in Appendix C. The underlying features of the fit in which the first two generations including the 
neutrinos are elementary from the ADS/CFT point of view while the third generation fermions 
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{tLiiR) having a tendency to be partially composite or composite, is maintained for both the SM 
and the SUSY case. From the choice of the c parameters, we find that the LLHH case admits a 
better fit to the neutrino data than the Dirac case. This is contrary to the observations made in 
[30] in normal RS where the c parameters for all the leptons were close to unity. It thus offered a 
more viable alternative than the LLHH case. In order to compensate for the weak warp factor in 
the Dirac case the right handed neutrinos had bulk masses cn ~ 7. This weak warping however, 
works in favour of the LLHH case where for c < 0.5 the effective 4D suppression scale is of the 
0{Mgut) resulting in fits with c parameters closer to unity. 

IV. SUSY SPECTRUM AND FLAVOUR PHENOMENOLOGY 

There are several ways to break supersymmetry within this RS set up at the GUT scale (see 
for example discussion in [31, 34-36, 52]. In the present work, we will consider only one particular 
set up which manifestly demonstrates the flavour structure of the fermions within the soft terms. 
More detailed analysis of supersymmetric spectrum will be addressed in [53]. We will assume in 
the following that supersymmetric breaking happens on the IR brane, or the GUT brane. Unlike 
the work of [33] and [36] we will not arrange the SM fields in any particular GUT representation. 
As has been discussed in these works, a GUT structure can be arranged with possible solutions for 
proton decay and doublet triplet splitting. Instead we parameterize SUSY breaking in terms of a 
single four dimensional spurion chiral superfield, X = 9'^F, which is localized on the GUT brane. 
However, for soft masses generated at the Planck brane as in [31], one may then impose the GS 
anomaly cancellation conditions on bulk masses to ensure unification of couplings at the Planck 
scale. We do not impose any such conditions as the soft masses are generated at the GUT scale. 

In the limit, the higher KK modes are decoupled from the GUT scale physics [38], the Kahler 
potential relevant for the scalar mass terms is given by 

/C(^) = / dy6{y-7rR)e-^''^'^k-^X^X [p,,^jQlQj + l3u,ijUJUj + |3d,^JDlDj + -fu,dHlaHu,d + ■■) 

(15) 

where /3 have dimensional carrying negative mass dimensions of -1 ( as the matter fields are five 
dimensional). ju,d s^-ie 0(1) parameters. 

The sfermion mass matrix is generated when the X fields get a vacuum expectation value 



m| r^ k ^ < X >' < X > Q'Q. The mass matrix will however not be diagonal in flavour space. 
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In the canonical basis, (15), the mass matrices take the form 

{^}h = ^1/2 k e(i-^-^^)^^-e(c.)e(c,) (16) 

where /3jj = 2kf3ij are dimensionless 0{1) parameters, ^(cj) are defined in Eq.(4). And the gravitino 
mass is defined as 



"^3/2 = -^^ = -Jar (17) 



The Higgs fields are localised on the GUT brane, their masses are given by m\j ^ = ^u4 "T'3/2- 

The A-terms are generated from the higher dimensional operators in the super potential of the 
type : 

VF(^) = j dy5{y - T:R)e-^^yk-^X [A^t^HuQ^Uj + Af^HaQidj + A^^HaLiEj + . . .) (18) 

where the A are dimensionful parameters having mass dimension -1. Substituting for the vev of 
the X, we have for the four dimensional trilinear couplings at the GUT scale: 

A^f = my,A[/'-^^-^'.^'''-ac^W,) (19) 

where we defined the dimensionless 0(1) parameters as A' = 2kA. The structure of the A terms 
and the corresponding fermion mass matrix are similar and they differ only by the choice of the 
0(1) parameters. Choosing A' = 2kY' , makes the down sector A terms diagonal in the mass basis 
of the fermions at the GUT scale. Henceforth, we shall work in this basis, with the 0(1) parameters 
of the A terms proportional to the 0(1) Yukawa parameters. 

The masses for the gauginos are obtained from the following operator in the lagrangian 

C= f d^ek'^XWAa'^A (20) 

At Mgut their masses will be be w-x /2 = fm^j2 where / is a 0(1) parameter, m^ /2 will be treated 
as an independent parameter. They are independent of the position of localization of X. as the 
profile for the gauginos is fiat corresponding to a bulk mass parameter of 0.5 [38, 40]. 

While the above equations set the boundary conditions at the high scale, the weak scale spec- 
trum is determined by the RGE evolution. In the present case, the spectrum at the high scale 
is completely non-universal as determined by the profiles of the zero modes of the matter chiral 
superfields. The structure of soft terms discussed here is similar to the ideas of flavourful super- 
symmetry discussed by [54, 55] and more recently by [56]. In the following we will present two 
example points one for the LHLH higher dimensional operator case and another for the Dirac case. 
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To begin with, in both the examples, we consider that all the 0(1) parameters appearing in the 
definitions of the soft parameters are proportional to the unit matrix. We will explicitly mention 
any deviations as required by the phenomenology when presenting numerical examples. This would 
mean that the matrices, A' , /3 in Eqs. (16, 19) are proportional to unit matrix and the parameters 
")u-,ld in Eq.(15) are equal to one. However, as we will see below, they play an important role 
in low energy phenomenology and one might frequently require to vary them within the 0(1) 
range, to satisfy phenomenological constraints. While studying the flavour phenomenology, we 
make sure that the soft terms are present in the super-CKM basis. The low-energy spectrum has 
been computed numerically using the spectrum generator SUSEFLAV [57]. 

3. LHLH operator case 

In this case we consider the following point, (21), in the c parameter space. It has a x^ of 5.5 
for the hadronic sector and 0.7341 for the leptonic sector. As expected it has a mostly composite 
right handed top quark. In addition, the leptonic doublets are also significantly composite in this 
case. 



cq^ = 2.740 CD, = 0.722 cu^ = 0.4024 cl, = -1.497 ce^ = 3.634 
CQ^ = 1.920 CD2 = 0.729 cu^ = 0.0652 cl^ = -0.224 ce^ = 2.290 
CQ3 = 0.960 CD3 = 0.801 ct/g = -3.5615 cl^ = -1.0738 ce,, = 1.769 (21) 

The choice of 0(1) parameters in the soft sector plays a role in determining the nature of the low 
energy spectrum. For a given set of c parameters, a naive choice of one for all the 0(1) parameters 
in the soft sector may or may not lead to an acceptable spectrum at Mgusy For the LLHH case 
corresponding to the choice in Eq.(21) the 0(1) parameters for all the soft masses are taken to 
be 1. The 0(1) parameters for the A terms are chosen to be A^ = 1.02y " while A" = Y and 
A^ = 0.6y ^. Corresponding to these choices of the 0(1) parameters and the c values in Eq.(21), 
the soft breaking terms at the GUT scale in GeV are given as: 



mq 



At, 



0.001 -0.03 -0.27 
-0.03 0.85 7.6 
-0.27 7.6 68.7 

-0.002 -0.09 -0.84 

-0.0002 1.09 -6.3 

-10-*^ 0.01 439.4 



\mu 



■.A 



D 



11.5 -63.5 156.1 

-63.5 349.6 -859.1 

156.1 -859.1 2110.8 

0.03 

-0.40 

-40.6 



;to_d = 



;mL 



105.03 -90.4 155.2 
-90.4 77.8 -133.6 
155.2 -133.6 229.5 

7.7 8.8 -147.5 

8.8 9.9 -167.0 
-147.5 -167.0 2798.4 
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ruE 



0.001 -0.017 0.05 

-0.017 0.26 -0.89 

0.05 -0.89 3.04 



]Ae = 



0.008 
1.81 
-31.3 



(22) 



A couple of interesting features of the above spectrum are (i) at least one of the soft masses 
is tachyonic (ii) significant amount of flavour violation present at the high scale. However at the 
weak scale, things are significantly different. This is because the RG running is quite different for 
the diagonal terms compared to the off-diagonal ones. In fact, the off-diagonal entries barely run, 
where as the corrections to the diagonal ones are quite significant. As an illustration, consider the 
slepton mass matrix at the weak scale. The analytic output at the weak scale for the diagonal 
terms can be approximated as 






(23) 



mi^{m 



susyj 



which receive gauge contributions while the off diagonal elements do not. The A terms are not 
large enough to make the off diagonal elements of the soft mass matrices comparable with the 
diagonal terms. An example for the sleptons for the case under consideration is given as 

2.2 X 10^ 77.0 -2.1 X 10^ 

77.0 2.2 X 10^ -2.7 x 10^ GeV^ 

-2.1 X 10"^ -2.7 X 10"^ 7.7 x 10^ 

1.1 X 10^ 3.4 X 10"^ 2.2 X 10^1 
susyj - 3.4 X 10""^ 1.1 X 10^ 5.9 X 10^ GeV^ 
2.2 X 10-1 5.9 X 10^ 5.4 x 10^ 

We see that the off-diagonal entry has barely enhanced where as the diagonal entries have been 
significantly modified. Constraints from fiavour violation would restrict the mass scales of m^/2 and 
Ml/2- The most stringent constraints are from the transitions between the first two generations ie 
from K^ —7- K^ and /i — )■ e-l-7. The expressions for the Kl — Kg mass difference and the branching 
fractions for /i — )■ e + 7 can be found in [58, 59]. We impose the fiavour constraints from all the 
existing data on the 5 parameters. 



m^^m. 



(24) 
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TABLE XI: Experimental upper bounds on the (5''°"'" obtained for rh„ = 2.1 TeV and rhi = 0.7 TeV 



(iJ) 


\s?,\ 


\sL\ 


\SEr\ 


l-^fi^l 


Kl\ 


I-^IlI 


\sRr\ 


Kr\ 


12 


0.053 


0.0002 


0.0003 


3.8x10-6 


0.0003 


3.8x10-6 


0.03 


0.03 


13 


0.34 


0.14 


0.06 


0.03 


0.06 


0.03 


0.26 


- 


23 


0.61 


0.16 


0.01 


0.04 


0.02 


0.04 


0.84 


- 



Bounds from the flavour violating processes are obtained using the mass insertion approximation 
[58] and the results of [60] defined in the Super-CKM basis. The flavour violating indices are defined 

as Sij{i 7^ j) = ^^ — —{i 7^ j) are evaluated in the basis in which the down sector is diagonal, 

with U being the rotation matrix which rotates the corresponding fermion mass matrix. The S 
are evaluated at the weak scale and we scale the bounds of [60] to the present mass scales. In 
Table[XII] we present the low energy spectrum corresponding to the sample point in Eq.(21). The 
low energy 6's are presented in Table XIII. 

TABLE XII: Soft spectrum for LLHH case: nisusy = 1-06 TeV, rrig = 2.64 TeV, ^i = 3.43TeV, tan/S = 25 



Parameter 


Mass(TeV) 


Parameter 


Mass (TeV) 


Parameter 


Mass (TeV) 


Parameter 


Mass(Tev) 


Parameter 


Mass (TeV) 


ii 


0.47 


6i 


1.01 


Tl 


0.726 


i>T 


2.78 


^1 


0.465 


i2 


1.05 


&2 


2.14 


T2 


2.79 


/>M 


0.483 


^2 


0.929 


cr 


2.24 


SR 


2.40 


fj-R 


0.478 


Ve 


0.469 


A^3 


3.38 


cl 


2.48 


SL 


2.48 


ML 


1.05 


- 


- 


A^4 


3.39 


ur 


2.24 


Ar, 


2.40 


SR 


0.476 


- 


- 


Cl 


0.895 


ul 


2.48 


dL 


2.48 


sl 


1.05 


- 


- 


C2 


3.43 


ruAQ 


3.23 


m± 


3.23 


ruh 


0.12186 


rriH 


3.06 


- 


- 



TABLE XIII: Low energy 5' s for quarks and leptons corresponding to the points in Eq.(21) for the LLHH 
case evaluated for fhq ^ 2.1TeV and m; = 0.7 TeV 



(iJ) 


l'^?.! 


I^^lI 


\SEr\ 


Kr\ 


\SRl\ 


Kl\ 


Kr\ 


Kr\ 


Kr\ 


12 


0.0003 


0.0001 


10-1" 


10-8 


10-8 


10-5 


0.001 


10-10 


0.001 


13 


0.01 


0.04 


io-« 


io-« 


10-5 


0.01 


0.005 


10-^ 


0.01 


23 


0.05 


0.05 


10-6 


10-5 


10-4 


0.04 


0.003 


0.0001 


0.07 
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4- Dirac Case 

For the case where neutrinos are of Dirac type the c parameters in Eq.(25) with x^ of 0.3211 for 
the hadronic sector and 0.1481 for the leptonic sector were chosen. The c values for the doublets 
in this case indicate they are predominantly elementary from the CFT point of view especially for 
the first two generations. The third generation however may be partially composite as in this case. 

CQ^ = 1.895 CD, = 1.898 cu, = 1.738 cl, = 1.293 ce, = 2.480 ctv, = 6.783 
CQ^ = 1.467 CD2 = 1-271 CU2 = 1-124 cl^ = 1-311 ce^ = 1-406 cn^ = 7-346 
CQ3 = -0.137 CDs = 1-394 cus = -0.356 c^g = 0.260 ce., = 0.237 cn, = 7.332 (25) 

The generic feature of soft mass matrices, discussed in the LLHH case, of the spectrum being 
tachyonic at the high scale and the diagonal terms evolving more than the off diagonal elements 
apply to this case as well. Corresponding to the c values in Eq.(25), the soft masses at the GUT 
scale have been evaluated for rn^V^ = 800 GeV while choosing Mi^ = 1200 GeV for the three 
gauginos. The 0(1) parameters corresponding to (w,q)33 and (m[/)33 were chosen to be 4 while for 
the others they are set to be 1. A^ = 1.151^ '^\/i,j, for all the A terms of the up sector while for the 
down sector and the leptons they were set equal to the corresponding 0(1) Yukawa couplings. The 
high scale soft breaking matrices in GeV are given in Eq.(26). In Table[XIV] we present the low 
energy spectrum corresponding to the sample point in Eq.(25). The low energy 5' s are presented 
in Table XV. 



mq 



0.60 1.2 30.0 
1.2 8.7 -19.4 
30.0 -19.4 2560.6 



■,mu 



0.7 -3.8 17.1 
-3.8 32.1 70.3 
17.1 70.3 2971.6 



■,mD 



0.10 -0.85 -1.7 
-0.85 7.1 14.6 
-1.7 14.6 29.8 



A 



u 



lO"-'^ -0.42 -5.8 
lO"-'^ 1.18 0.20 
10-5 -0.005 488.8 



:A 



D 



0.03 

0.57 
-40.6 



■,rnL 



8.7 0.9 62.5 
0.9 0.1 7.0 
62.5 7.0 446.3 



niE 



0.2 -0.4 10.5 

-0.4 0.7 -17.9 

1.0 -17.29 442.7 



■,Ae 



-0.01 

-3.01 
52.1 



(26) 
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TABLE XIV: Soft spectrum for Dirac case: rususy = 1-05 TeV, rrig = 2.65 TeV, ^ = 4.32TeV, tan[3 = 25 



Parameter 


Mass(TeV) 


Parameter 


Mass(TeV) 


Parameter 


Mass(TeV) 


Parameter 


Mass(Tev) 


Parameter 


Mass (TeV) 


ii 


0.702 


6i 


2.06 


Tl 


0.480 


Dr 


0.570 


^1 


0.465 


h 


2.31 


h 


2.32 


T2 


0.802 


/>M 


0.624 


N^ 


0.928 


cr 


2.25 


SR 


2.36 


fj-R 


0.608 


Ve 


0.625 


N3 


4.26 


cl 


2.45 


SL 


2.45 


Al 


0.902 


- 


- 


7V4 


4.26 


UR 


2.25 


Ar 


2.36 


SR 


0.610 


- 


- 


Cl 


0.894 


ul 


2.45 


dL 


2.45 


Cl 


0.903 


- 


- 


C2 


4.32 


niAQ 


4.18 


m± 


4.18 


mh 


0.1235 


niH 


3.96 


- 


- 



TABLE XV: Low energy S's for the Dirac Case corresponding to the point in Eq.(25) for Dirac case evaluated 
for rhq = 2.1TcV and m; = 0.7 TeV 



(ij) 


\s2,\ 


\sL\ 


\SSr\ 


\sU 


Kl\ 


Kl\ 


\5Rr\ 


1-51^1 


Kr\ 


12 


0.0003 


10^6 


10-10 


10-*^ 


10-s 


10-5 


10-7 


10-7 


0.00005 


13 


0.01 


0.007 


io-« 


10-s 


10-5 


0.002 


10-6 


10-4 


0.06 


23 


0.06 


10-4 


10-6 


10-5 


10-5 


0.01 


10-4 


0.0006 


0.001 



V. OUTLOOK 

The Randall-Sundrum framework is typically considered to be the geometric avatar of the 
Proggatt-Nielsen models. In the present work, we have considered a warped extra dimension close 
to the GUT scale. We fit the quark masses and the CKM mixing angles and determined the range 
of the c parameters which give a reasonable x^ fit. The 0(1) parameters associated with the 
Yukawa couplings have also been varied accordingly. Though the top quark Yukawa is smaller at 
the high scale compared to the weak scale, it is still large enough that one requires a large negative 
bulk mass parameter for the right handed top quark. For the leptons, we considered two particular 
models for neutrino masses (a) with Planck scale lepton number violating operator and (b) Dirac 
neutrino masses. 

The results show that there is a significant difference in the RS models at the weak scale and 
the RS models at the GUT scale especially if one focuses on the neutrino sector. In the weak scale 
models, the Planck scale lepton number violating higher dimensional operator was very hard to 



22 



accommodate with perturbative Yukawa couplings and thus was highly disfavoured. The Dirac and 
the Majorana cases were favoured though they were strongly constrained by the data from flavour 
violating rare decays. The situation is sort of reversed in the GUT scale RS models, though not 
exactly. The higher dimensional Planck scale operator fits the data very well, where as the Dirac 
case requires larger c values, some times with bulk mass parameters almost an order of magnitude 
larger than the cut-off scale. In the hadronic sector, the situation is not so dramatic. The top 
quark Yukawa though reduces at the GUT scale compared to the weak scale, still requires that the 
right handed top to be located close to the IR, making it a composite as in the weak scale models. 
The supersymmetric version of the same set up is far more interesting as it can have possible 
observable signatures at the weak scale. The main difference in fitting of the fermion masses is due 
to the presence of the additional parameter tan/? in the supersymmetric case. The difference in the 
c values is not very significant for low values of tan (3. At large tan /3, for a given generation, the zero 
modes are localized more towards the IR brane as compared to the SM case. This effect is more 
pronounced in the down sector as shown in Figures [7]. We parameterise SUSY breaking by a single 
spurion field localised on the IR brane. The resultant soft masses depend on the profiles of the zero 
modes of the chiral superfields and contain fiavour violation. At the weak scale the constraints 
from first two generation fiavour transitions rule out light spectrum. Another significant constraint 
comes the light Higgs mass at 125 GeV. The trilinear couplings have the same form as the Yukawa 
couplings in this model as long as the 0(1) parameters associated with both the parameters are 
taken to be proportional to each other. If all the 0(1) parameters at high scale are take to be 
exactly unity, the weak scale values of At are small to generate a 125 GeV light Higgs, for stops 
of masses ~ 1.5 — 2 TeV. However, with a minor variation of the 0(1) parameters for At at the 
high scale, 125 GeV Higgs is easily possible. More variations of supersymmetric breaking and the 
corresponding spectra will be discussed later [53]. 
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APPENDIX A: PLOTS FOR RANGES OF C PARAMETERS FOR QUARKS AND LEP- 
TONS FOR SM FITS AT THE GUT SCALE 

1. Range of c parameters for the quarks 
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FIG. 1: The points in the above figures correspond to a x^ between 1 and 10 for the SM fits at the GUT 
scale. The plot represents the parameter space for the bulk masses of the quarks. 
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2. Range of c parameters for the for leptons for the LLHH scenario 
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FIG. 2: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the LLHH case. 
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3. Range of c parameters for the leptons. 
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FIG. 3: The points in the above figures correspond to a x^ between 1 and 10 for the SM fits 
scale. The plot represents the parameter space for the bulk masses of the leptons corresponding 
type neutrinos. 
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APPENDIX B: PLOTS FOR RANGES OF C PARAMETERS FOR QUARKS AND LEP- 
TONS FOR THE SUPERSYMMETRIC CASE 

1. Range of c parameters for the quarks 
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FIG. 4: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the quarks. 
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2. Range of c parameters for the for leptons for the LLHH scenario 
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FIG. 5: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the LLHH case. 
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3. Range of c parameters for the for leptons for the case of Dirac neutrinos. 
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FIG. 6: The points in the above figures correspond to a x^ between 1 and 10. The plot represents the 
parameter space for the bulk masses of the leptons corresponding to the Dirac type neutrinos. 
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APPENDIX C: COMPARATIVE PLOT BETWEEN SM AND SUPERSYMMETRIC 
FITS FOR FERMIONS 
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FIG. 7: The plots correspond to eomparison between SM and supersymmetric fits for down sector fermions. 
Red corresponds to SM, wliile Blue and Green corresponds to supersymmetric case for tan/3 = 10 and 
tan/3 — 50 respectively 
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